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ABSTRACT 

We present the results of near-infrared multi-object spectroscopic observations for 37 Bz/ir-color-selected 
star-forming galaxies conducted with MOIRCS on the Subaru Telescope. The sample is drawn from the Kg- 
band selected catalog of the MOIRCS Deep Survey (MODS) in the GOODS-N region. About half of our 
samples are selected from the publicly available 24 /im-source catalog of the Multiband Imaging Photometer 
for Spitzer on board the Spitzer Space Telescope. Ha emission lines are detected from 23 galaxies, of which 
the median redshift is 2.12. We derived the star formation rates (SFRs) from extinction-corrected Ha lumi- 
nosities. The extinction correction is estimated from the SED fitting of multi-band photometric data covering 
UV to near-infrared wavelengths. The Balmer decrement of the stacked emission lines shows that the amount 
of extinction for the ionized gas is larger than that for the stellar continuum. From a comparison of the extinc- 
tion corrected Ha luminosity and other SFR indicators we found that the relation between the dust properties 
of stellar continuum and ionized gas is different depending on the intrinsic SFR (differential extinction). We 
compared SFRs estimated from extinction corrected Ha luminosities with stellar masses estimated from SED 
fitting. The comparison shows no correlation between SFR and stellar mass. Some galaxies with stellar mass 
smaller than -- lO'^M© show SFRs higher than -- lOOMoyr"^. The specific SFRs (SSFRs) of these galaxies 
are remarkably high; galaxies which have SSFR higher than ^ 10"^ yr"' are found in 8 of the present sample. 
From the best-fit parameters of SED fitting for these high SSFR galaxies, we find that the average age of the stel- 
lar population is younger than 100 Myr, which is consistent with the implied high SSFR. The large SFR implies 
the possibility that the high SSFR galaxies significantly contribute to the cosmic SFR density of the universe at 
z^2. When we apply the larger extinction correction for the ionized gas or the differential extinction correc- 
tion, the total SFR density estimated from the Ha emission line galaxies is 0.089 -0.1 36 MQyr"'Mpc~^, which 
is consistent with the total SFR densities in the literature. The metallicity of the high-SSFR galaxies, which 
is estimated from th e N2 index, is lar ger than that expected from the mass-metallicity relation of UV-selected 
galaxies at z ^ 2 by Erb et al.l ( l2006ah . 
Subject headings: galaxies: evolution — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

Recent studies revealed that a significant fraction of stars 
in the present-day galaxies were formed between redshifts 
1 and 3 (e.g., j Hopkins 2004; Perez-Gonzalez et al. 2005; 
Hopkins & Beacom 2006; ,Wang et al, .2006; .Dahlen et al, 
2007; Caputi et al. 2007). For example, results of deep near- 
infrared observations show the average stellar mass density 
at 0.5 < z < 1 is ^ 53-72% of that in the present universe, 
while at 2.5 < z < 3.5 it is just ^ 4-9% of the present value 
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jKaiisawa et al.l 12009^. Although such studies demonstrate 
the total star formation rate (SFR) density in the universe as 
a function of cosmic time, it is still uncertain which popu- 
lations of galaxies contribute to the active star formation be- 
tween redshift 1 and 3. The evolution of the stellar mass-SFR 
relation constrains theoretical views of how galaxies accumu- 
late the stellar mass (Dave 2008). Therefore, as a next step, in 
order to understand the build up of stellar mass in galaxies at 
this redshift range, it is crucial to examine the relationship of 
SFR to stellar mass. 

SFRs are often estimated from photometric indicators 
that are closely related to the emission from massive stars 
jKennic utt 1998). Massive stars have short lifetimes, so 
they represent stars which were born recently. The UV con- 
tinuum (1500-2800 A) directly reflects the number of mas- 
sive stars, and therefore can be used as an indicator of SFR 
for galaxies in this redshift range. Additionally, far-infrared 
emission is another indicator of SFR; the UV emission of 
young massive stars is absorbed by interstellar dust and rera- 
diated in the far-infrared at ^8-1000 /im. Fluxes at 24 /im 
wavelength taken with the Multiband Imaging Photometer for 
Spitzer (MIPS) are commonly used to estimate far-infrared 
emission of galaxies at high redshifts, because observations 
in lo nger wavelengths are too shallow to detect such galax- 
ies. iDaddi et"al] (l2007ah showed that /T-selected galaxies 
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at z ~ 2 have a tight correlation (^ 0.2 dex) between stel- 
lar mass and SFR estimated from UV and IR luminosities. 
The SFRs of these galaxies are roughly proportional to their 
stellar masses. However, the relation is not still conclusive. 
Other studies show different relationships between SFRs and 
stellar masses. For example, rest-UV selected BM/BX galax- 
ies also have a correlation bet ween these values b ut the cor- 
relation is not proportional tReddy et al.| l2006ah . Further- 
more, the SFRs of 24 um-selected galaxies do not depend on 
their stellar masses (iCaputi et al.ll2006a ). We need to note 
that SFRs estimated from UV and IR luminosities have large 
uncertainties. UV light is affected by dust extinction and 
stars with long lifetimes of ^ 1 Gyr. On the other hand, 
24 fim radiation corresponds to the rest-frame 8 /im wave- 
length for galaxies at z ^^ 2. It is known that the far-infrared 
luminosity estimated from rest-8 /im luminosity has a large 
uncertainty due to the complex SED at those wavelengths 
JPeeters et al. 2004; Chary & Elbaz 2001: Daddi et al. 20071 
[Papovich et al. 2007). 

Ha luminosity is one of the most reliable indicators of SFR. 
First, the emission line is reemitted from hydrogen gas ion- 
ized by UV photons shorter than 912 A, which are radiated 
by only the most massive stars. These stars have shorter life- 
times than stars emitting the UV continuum above 1500 A. 
Second, the emission line is less affected by dust extinction 
than UV emission. The amount of dust extinction is larger at 
shorter wavelengths. Third, the Ha luminosity can be mea- 
sured without the uncertainty inherent in SED estimates for 
the far-infrared luminosity. However, it is difficult to mea- 
sure Ha luminosities of galaxies at high redshift, because the 
Ha emission line at z > 0.4 is out of the wavelength range of 
optical spectrographs. 

Recently, near-infrared spectroscopic studies of SFRs based 
on the Ha luminosity for large sampl es of star- forming galax- 
ies at z ~ 2 have b een conducted JErb et al .l l2006bl . 120031: 
iHavashi et al]l2009t lF5rsFer Schreiber et al.l l2009). Er b et all 
(l2006bl) demonstrate a weak correlation between SFR esti- 
mated from Ha luminosity and stellar mass. Their sample 
galaxies are selected by U„GTZ color and have redshifts con- 
firmed using a blue-sensitive optical spectrograph. The ob- 
served galaxies are biased toward less attenuated galaxies. 
Since red galaxies account for the s ignificant fraction of stel- 
lar mass density at z '^ 2 (iPerez-G onzalez et al. 2008), the 
galaxies at z '^ 2 should be re-examined with a technique 
more sensitive to red galaxies (e.g., /T-selected galaxies). The 
A"-sele cted star-forming galaxies observed bv IHavashi et al.l 
( 120091) showed a flatte r SFR-stellar mass correlation than that 
of lDaddi et al.l (l2007ab . but the number of /T-selected galaxies 
with Ha detection is still rather limited and observations for a 
larger sample are necessary. 

In this work, we study the SFR and stellar mass of a large 
number of /T-selected star-forming galaxies at z ~ 2 based 
on near-infrared observations with the Multi-Obi ect InfraRe d 
Camera and Specti'ograph (MOIRCS; ISuzuki et al.1 l2008h . 
MOIRCS has the capability for multi-object near-infrared 
spectroscopy. We also utilize the deep near-infrared imaging 
of the MOIRCS Deep Survey (MODS; Kajisawa et al. 2009), 
which covers ~ 103.3 arcmin" in the Great Observatories Ori- 
gins Deep Survey North (GOODS-N) region. The sample for 
the present work is selected from ^-selected galaxies in the 
MODS catalog and their properties are studied using MODS 
and publicly available multi-wavelength data of GOODS-N. 
Throughout the paper, we use the AB magnitude system and 



adopt cosmological parameters of h = //q/IOO = 0.7, il^f = 
0.3, f7A = 0.7. 

2. OBSERVATION 
2.1. Sample Selection 

In order to select star-forming galaxies at z ~ 2, we ap- 
ply the BzK color criterion for star- forming galaxies (sBzK; 
BzK = {z-K)-(B-z) > -0.2; iDaddietal.) 12004) to the K- 
selected galaxies with Ks < 23 from the MODS catalog, 
whose limiting magnitude is Ks = 24.9 (5cr, (/>1.2arcsec). The 
BzK diagram of the extracted galaxies is shown in Figure [T] 
B and z band magnitudes are measured on the F435W and 
F850LP images of the Advanced Camera for Surveys (ACS) 
on board the Hubble Space Telescope (HST) in the course of 
the G OODS Treasury Program (Version 1.0; lGiavalisco et akl 
I2OOJ). In total, 508 galaxies are identified as candidates. 

We chose 37 sBzK galaxies from the candidates for the 
spectroscopic observation. The coordinates and photometry 
of the target galaxies are listed in Table [T] Although two 
of the sample galaxies (MODS 12-0125 and MODS41-0194) 
are not detected in the B-band, we include these galaxies be- 
cause their z-K colors suggest that they are sBzK or passive 
BzK (pBzK; BzK > -0.2 f]z-K > 2.5), and their photometric 
redshifts are consistent with the redshift range of BzK galax- 
ies. We did not select on the basis of spectroscopic redshift, 
although spectroscopic redshifts were previously known for 
some galaxies in the present sample (TablefTll. 

To preferentially select actively star-forming galaxies, we 
put higher priority on the galaxies identified in the publicly 
available Spitzer MIPS source catalog (the first delivery of the 
Spitzer data catalog from the GOODS Legacy Project; Chary 
et al. in preparation). The flux limit of the catalog (524 > 
80 yitJy) corresponds to L7-/R-- 4.1 x 10^'-3.6 x W^hp, atz = 
1.4-2 .5 using the infrared SED templates bv lCharv & Elbazl 
(|2001"); thus the galaxies in the MIPS source catalog are as 
luminous as (ultra-)luminous infrared galaxies. MIPS sources 
in the public catalog are identified with 127 of the sample 
galaxies (we call them sBzK-MIPS galaxies, hereafter). We 
selected 18 sBzK-MIPS galaxies. Furthermore, we added 19 
sBzK galaxies not listed in the MIPS public catalog to the list 
of targets (sBzK-non-MIPS galaxies, hereafter). 

It is known that sBzK criterion picks out various types 
of star-formin g galaxies and pass ively evolving galaxies at 
1 .4 < z < 2.5 fcrazian et alJl2007h . In order to check that our 
sample represents the sBzK population, we compare the B-K 
color distribution of our target galaxies to that of sBzK-non- 
MIPS and sBzK-MIPS galaxies with K^ < 23 in the MODS 
field (Figure |2]i. The B-K color reflects both the age and 
the amount of dust of a galaxy at this redshift range. Our 
sBzK-MIPS sample is reasonable representative of those in 
the field, according to both x^-test (64.65 %) and Kolmorov- 
Smirnov (KS) test (92.46 %). On the other hand, the sBzK- 
non-MIPS sample is possibly biased toward blue galaxies 
(X^-test : 35.35 %, KS-test: 17.74%). Red sBzK-non-MIPS 
galaxies are relatively rare, so that the scaled number in the 
redder bin is too small (< 1) to provide robust statistics given 
the present sample size. 

We utilize publicly available multi-wavelength data of 
GOODS-N to exclude possible AGN which may contami- 
nate the optical star-formation emission. The hard X-ray lu- 
minosity at the detection limit of the Chandra Deep Fie ld- 
North 2 Ms Point-Source Catalogs (lAlexander et al.ll2003h at 
z ~ 2 (1 X lO'^^ergs"'), suggests that the emission lines of 
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Fig. 1. — BzKcolor-colordiagramof the a;, < 23 objects from the A^^-band selected catalog of the MOIRCS Deep Survey (MODS: lKausawa et aU2009l) . Red 
syinbols show the MIPS sources listed in the pubHc catalog (Chary et al. in preparation). The sample for the present work is denoted by circles and triangles for 
the objects with emission lines detected and not detected, respectively. Arrows show upper/lower limits for sBzK galaxies. 
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Fig. 2. — B—K color distribution of sBzK-MIPS (left) and sBzK-non-MIPS (right) galaxies (see text). The bar plots show our observed sample. The dashed 
line with error bars shows the scaled color distribution of all of the Ks < 23 sBzK-MIPS and sBzK-non-MIPS galaxies, respectively. The error bars show the 
Poisson eiTors at each bins. 
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the hard X-ray detected sBzK galaxies are powered by an 
AGN rather than massive stars ( [Alexander et al.ll2002b . X-ray 
sources are identified with 62 of the 508 sBzK galaxies. We 
exclude these sources from the present sample. Rest-frame 
infrared colors are often used to identify AGN candidates us- 
ing Spitzer IRAC bands (3 .6 /im, 4.5 iim, 5.8 M m, and 8.0 fim) 
(e.g.. lCaputi"era l. 2006b; Havashi et al.ll2009l) . A galaxy with 
emission dominated by the stellar population has a strong 
bump in the SED at 1.6 /im JSawi cki 2002). We confirmed 
that all of the present non-X-ray galaxies have strong 1.6 //m 
bumps in the SED of the IRAC bands. 

2.2. Observations 

The spectroscopic observatio ns for the selected galaxies 
were performed with M OIRCS ([Suzuki et aI1l2008 ') attached 
to the Subai-u Telescope C ive et alj|2004l) during 2007 Mai'ch 
24-27 in the course of the near-infrared spectroscopic obser- 
vation of X-ray sources (S07A-083, PI M. Akiyama). We 
used MOIRCS in MOS mode with the low-resolution grism 
(HK500), which covers 1.3-2.3 /im wavelengths. The spec- 
tral resolution is /? ~ 500 (Av -- 600 km s"') with 0''8 slit 
width. The 37 galaxies were observed with four slit masks 
as shown in Table |2] One slit mask was used for each night. 
The total exposure time was 2.5-5 hour for each mask. The 
slit width was fixed to 0"8 for targets. Most of the slit 
lengths were 10. "0 or longer, depending on target positions. 
To monitor the flux variation due to the atmospheric absorp- 
tion and photon loss at the slit mask, a bright point source 
was observed simultaneously on each mask, except for mask 
CDFN3. 

The FWHMs of seeing size in the K, band during the obser- 
vations were typically ^ 0"6. According to the atmospheric 
attenuation archived in the Canada France Hawaii Telescope 
Web siteH, the conditions were photometric for the first three 
nights. For the forth night, thin cirrus covered the sky. The 
observation log is summarized in Table |2] 

The observation followed the standard proced ure for MOS 
observation with MOIRCS described in Tokokul (I2007h . The 
target acquisition was done using 6 to 7 alignment stars dis- 
tributed over each slit mask. The telescope was dithered along 
the slit (we call it an 'AB pattern") by 3."0, to achieve suffi- 
cient sky subtraction. The detector was read eight times in 
non-destructive readout mode to reduce the read out noise 
down to 16 e~. The exposure time was set so that most of 
sky emission lines, which were used later for the wavelength 
calibration, were not saturated. The typical exposure times 
were 300-900 sec. The relative position of the alignment star 
against the slit mask was checked, and then the telescope 
pointing was adjusted every hour or two without removing 
the slit masks from the focal plane. The systematic guide 
error in one hour was measured to be <0."4. At the begin- 
ning of each observing night, an AOV type star (HIP53735; 
Kg = 8.856; Two Micron All Sky Survey point source catalog; 
ISkrutskie et al.ll2006.) was observed for flux calibration with 
the same instrument configuration as the targets. 

3. DATA REDUCTION AND DATA ANALYSIS 

3.1. Data Reduction 

The data reduction was made, using a semi-automatic 
reduction script (MOIRCS MOS Data Pipelines; MC- 
SMDlO)> which follows the standard procedure of 



IRAlQONEDSPEC. At first, we make a map of bad 
pixels; bad pixels are identified as those that show non- 
linearity. This includes dead (always zero count) and hot 
(always high count) pixels. Cosmic rays are detected by a 
moving block average filter using IRAF:CRAVERAGE. Each 
input frame for IRAF:CRAVERAGE is first divided by its 
dithered pair so that cosmic rays stand out behind the sky 
background. The cosmic -ray mask was combined with the 
bad pixel mask. The values of bad pixels and cosmic -ray 
pixels are replaced with values interpolated linearly from 
the pixels along the slit. Sky emission lines were subtracted 
using the dithered pair frames (we call the procedure 'A-B 
subtraction"). The dome-flat frames for flat fielding were 
taken with the same slit mask and grism configuration as used 
for the observation. For the first-order distortion correction, 
the database for the imaging mode of MOIRCS was applied 
to the spectral frames. The tilt of the spectrum against 
the detector rows was corrected in the cases when it was 
measurable in the continuum of bright objects. 

Next, each spectrum was cut out based on a mask de- 
sign file. For wavelength calibration, OH sky emission lines 
were identified in each frame using a night-sky spectral atlas 
(Rousselot et al. 2000). Using the pixel coordinates of the sky 
emission lines, the sky-subtracted frames were transformed 
such that the sky lines aligned along the column of the frames 
and that the wavelength was a linear function of pixel coordi- 
nates. Background emission still remained after A-B subtrac- 
tion due to time variation of OH emission during dithering 
observations. The residuals were subtracted by fitting with 
a quadratic equation along the column. Finally, all frames 
were weighted with the flux of the bright point source and 
stacked using the weighted average to obtain the final two- 
dimensional spectra. 

The response spectra, including atmospheric transmission 
and the instrument efficiency, are calculated using the flux 
standard star to convert the observed count to an absolute 
flux. The data of the flux standard star is reduced to a two- 
dimensional spectrum in the same way for galaxy targets and 
then the continuum is traced with a polynomial to extract a 
one-dimensional spectrum with IRAF:APALL. We calculate 
a model spectrum of the AO V flux standar d star using SPEC- 
TRUIVQ (SPECTRUM uses lCastelli &Kurucz 2004 model). 
The model was scaled to Ks = 8.856 and the spectral resolu- 
tion of the model was degraded to R ^ 500 with a Gaussian 
kernel, so as to give the same hydrogen line width as the ob- 
served spectra. Finally, the one-dimensional spectrum of the 
standard star was divided by the model spectrum to obtain the 
response spectrum. 

Emission lines are detected by visual inspection on the flux- 
calibrated two-dimensional spectra. The signal-to-noise ratio 
evaluation described later shows that visual inspection corre- 
sponds to a larger than ~ 3cr detection. The detected emission 
lines were integrated along slit direction within an appropri- 
ate aperture giving the largest S/N to obtain one dimensional 
spectra. 

3.2. Measurement of Emission Lines 

We first consider the strongest emission line as Ha. Taking 
account of the wavelength range where we observed and the 



'^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 

omy (AURA) under cooperative agreement with the National Science Foun- 

^'' ihttp : //nenue ■ cfht ■ hawaii ■ edu/Instruments/Elixir/sk yp Jafaai/home . html 

^' available at |http: //www.nao j ■ org/Observing/DataReduction/ '^ http: //ww w.phys ■ appstate ■ edu/spectrum/spectrum.html| 
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TABLE 1 
Properties of Galaxies Observed 



ID 


RA 


Dec 


K," 


B-z" 


z-K" 


S24' 


MASK ID 


^literature 


Zphot 


ZHa" 


emission line^ 




(J2000.0) 


(J2000.0) 


(mag) 


(mag) 


(mag) 


(MJy) 












MODS 11-0076 


12 36 18.50 


62 09 03.5 


20.5 


1.9 


2.5 


130.1 


CDFN2 


_ 


J 4g-tfl.02 

T 1 i+i);oi 

-0 18 


1.6731 


Ha 


MODS 11 -0094 


12 36 24.29 


62 08 27.9 


22.2 


1.0 


1.2 


106.8 


CDFN2 


- 


1.9976 


Ha, [Oin] 


MODS 11-0274 


12 36 25.29 


62 10 35.6 


21.4 


1.5 


2.4 


106.4 


CDFN2 


- 


2 ^n*t|'■^o 


2.0808 


Ha 


MODS 11-0299 


12 36 35.64 


62 09 41.4 


22.9 


1.2 


1.3 


38.6 


CDFN2 


- 


7 97+0:06 

-0.06 


- 


- 


MODS 11-0390 


12 36 39.37 


62 10 06.6 


21.8 


1.5 


1.5 


127.2 


CDFN2 


- 


2.3516 


Ha, H/3, [OIII] 


MODS 12-0 125 


12 36 11.50 


62 10 33.7 


21.7 


> 1.5 


3.1 


282.1 


CDFN2 


- 


2 j2"*06 
7 Ti+o:og 

"^■-^^-0 06 


2.2461 


Ha, [Nil] 


MODS 12-0255 


12 36 18.38 


62 11 39.1 


22.8 


1.1 


1.4 


9.7 


CDFN2 


2.390 


2.3979 


Ha, H/3, [OIII] 


MODS2 1-2005 


12 37 01.04 


62 10 16.8 


22.5 


1.0 


2.0 


210.7 


CDFN3 


- 


1 98+"-''' 


- 


- 


MODS21-2022 


12 36 51.80 


62 11 14.8 


22.6 


0.7 


0.8 


0.0 


CDFN3 


- 


- 


- 


MODS2 1-2047 
MODS21-2321 


12 36 54.95 
12 36 44.65 


62 10 56.4 
62 12 27.5 


22.0 
22.3 


1.4 
0.8 


2.0 

1.3 


204.2 
26.6 


CDFN3 
CDFN3 


0.999 
1.879 


2 j7+fl.05 

1 56+8:8§ 
1 cc+o:o7 


1.7341 


Ha 


MODS21-2612 


12 37 01.07 


62 10 54.3 


22.1 


1.9 


2.3 


99.0 


CDFN3 


- 


1.5256 


Ha 


MODS2 1-4647 
MODS21-5277 


12 37 11.85 
12 37 02.91 


62 11 15.6 
62 14 04.6 


22.4 
21.9 


0.7 

I.l 


0.9 

1.2 


7.2 
97.6 


CDFN3 
CDFN3 


1.244 


J 4j-tfl.06 

1 24+8:8i 
221+8:?? 

^■^^-0 52 


- 


- 


MODS21-6131 


12 37 08.77 


62 12 57.9 


22.6 


0.5 


1.0 


12.6 


CDFN3 


2.269 


2.2678 


Ha 


MODS22-1411 


12 36 26.97 


62 13 17.3 


22.0 


1.7 


2.1 


148.3 


CDFN2 


- 


2 i9+fl.08 
2 j3+o:oi 
1 ^^+o:o6 


2.3983 


Ha, [Nil], H/3 


MODS22-2282 


12 36 21.74 


62 14 52.9 


22.6 


0.7 


1.2 


12.7 


CDFNl 


2.200 


2.2087 


Ha 


MODS22-2658 


12 36 30.08 


62 14 27.9 


20.6 


2.1 


2.1 


181.1 


CDFNl 


1.525 


1.5246 


Ha 


MODS22-4194 
MODS22-4806 


12 36 33.67 
12 36 57.22 


62 15 33.0 
62 14 29.9 


22.8 
22.2 


1.0 
1.0 


1.0 

1.8 


13.4 
114.2 


CDFNl 
CDFNl 


2.488 
1.873 


239+0.04 
1.73+8-'' 

^■"^-0.06 


2.4862 


Ha 


MODS22-5133 


12 36 45.84 


62 14 47.0 


22.5 


1.5 


1.9 


48.7 


CDFNl 


- 


1.5255 


Ha 


MODS22-5375 


12 36 50.11 


62 14 01.1 


22.8 


0.5 


1.4 


69.3 


CDFNl 


2.231 


2 23+O.04 

2.36+8:^^ 
1 4<r+o:o2 

^ -^^-0 01 


- 


- 


MODS3 1-0033 


12 37 20.05 


62 12 22.8 


21.9 


1.2 


1.8 


113.3 


CDFN3 


2.458 


2.4605 


Ha, [Nil], H/3, [OIII] 


MODS3 1-0129 


12 37 26.43 


62 13 30.3 


21.1 


2.1 


2.1 


89.2 


CDFN3 


- 


- 


- 


MODS3 1-0199 
MODS32-0035 


12 37 19.39 
12 36 37.88 


62 16 21.1 
62 16 45.3 


22.5 
21.7 


0.7 

1.5 


0.9 
2.0 


13.8 
117.8 


CDFN4 
CDFNl 


1.567 


1 4g+fl.04 

J 45+8:05 
1 44+8:04 


1.5681 


Ha 


MODS32-0116 


12 36 50.26 


62 16 55.7 


21.9 


1.2 


1.6 


54.0 


CDFNl 


1.315 


1.4878 


Ha 


MODS32-0153 


12 36 53.65 


62 17 24.3 


21.9 


0.7 


1.6 


71.6 


CDFNl 


2.186 


2 02+«.l5 

183^8:15 
7 34+o:o2 

^■-^ -0.06 


2.1865 


Ha, [OIII] 


MODS32-0162 


12 36 51.33 


62 17 50.9 


22.4 


0.8 


2.0 


158.5 


CDFNl 


- 


- 


- 


MODS41-0194 


12 37 34.42 


62 18 12.1 


21.9 


> 1.1 


3.4 


118.9 


CDFN4 


2.241 


2.2425 


Ha, [OIII] 


MODS4 1-0297 
MODS42-0112 


12 37 24.15 
12 37 16.34 


62 16 11.7 
62 19 20.5 


21.9 
22.8 


1.1 
0.6 


1.6 
1.0 


123.6 
19.8 


CDFN4 
CDFN4 


1.689 


1 46+0.03 

2 jg+o:o5 
7 3f;+o:o5 


1.6855 
2.1231 


Ha 
Ha, [Om] 


MODS42-0135 


12 37 20.04 


62 19 23.2 


22.1 


1.1 


1.8 


78.1 


CDFN4 


2.295 


2.3049 


Ha, [Om] 


MODS42-0145 


12 37 24.11 


62 19 04.9 


22.1 


0.5 


0.9 


32.0 


CDFN4 


2.094 


205+fl.08 

1 45+8:45 


2.0952 


Ha 


MODS42-0162 


12 37 25.80 


62 19 14.2 


20.8 


1.6 


2.2 


61.4 


CDFN4 


- 


- 


- 


MODS42-0194 


12 37 33.50 


62 20 01.2 


22.0 


1.5 


1.6 


51.6 


CDFN4 


- 


1 36+0.07 
1 Qo+8:?? 

^■^^-0.17 


- 


- 


MODS42-0228 


12 37 22.28 


62 20 39.1 


21.9 


1.2 


1.6 


244.3 


CDFN4 


- 


- 


- 



^ Total magnitude, MAG_AUTO by SEXTRACTOR 

Color in 0I."5 apeiture centered at K^ coordinate 
^ 24 /im flux measured from the Spitzer 24 ^m data (Dickinson et al. in preparation; see j4.3t 
** Spectroscopic redshift in literature ICohe n 2001: Cohen et al.200a ; .Dawson et al.20 0I;. Wii-th et alJ200HCo 

l2006tHB^ra5- et aU200ai ^^ 

' Photometric redshift IKaiisawa et aU2009n 

Our spectroscopic redshift derived from Ha emission line 
8 Emission hnes detected with S/N> 3. [N II] and [O III] is [N 1I]A6583 and [O III] A5007, respectively. 



i et alJ2004ljTreu et alJ2003:IChaDman et alJ200aiReddv et alj 



TABLE 2 

Summary of the observation 



mask ID R.A. 

(J2000.0) 



Dec. 

(J2000.0) 



P.A." 
(degree) 



UT date 



target 
(number) 



exposure seeing" 
(min) (arcsec) 



CDFNl 


12 36 38.7 


62 15 58 





2007 Mar 25 


HK500 


10 


316 


0.52 


CDFN2 


12 36 22.3 


62 10 38 


249 


2007 Mar 24 


HK500 


8 


160 


0.64 


CDFN3 


12 37 03.2 


62 12 20 


338 


2007 Mar 26 


HK500 


10 


310 


0.68 


CDFN4 


12 37 21.3 


62 18 28 


228 


2007 Mar 27 


HK500 


9 


274 


0.69 



' Direction of the slits, from noith to east 

' Average FWHM sizes of PSF in Kg band during the observations 
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MODS22-22S2 



MODS22-265S 



MODS22-4194 



MODS22-5133 



MODS32-0116 




2.02 2.04 2.06 

Fig. 3 . — [N II]-Ho spectra of the emission-line galaxies. Shaded region shows the size of the error as a function of wavelength. The dashed line shows the 
spectrum of the best-fit model. The vertical dotted lines in each panel show, from left to right, the locations of [N II] A6548, Ho, [N II]A6583. 



MODS32-0153 



MODSll-0094 




MODSll-0390 MODS12-0255 

0.6- 



MODS22-1411 




gil.2 



1.54 1.56 1.58 1.60 t.44 1.46 1.48 1.50 

MODS31-0033 MODS41-0194 



1.62 1.64 1.66 1.68 1.64 1.66 1.68 1.70 1.64 1.66 1.68 1.70 

MODS42-0112 MODS42-0135 




1.68 1.70 1.72 1.74 1.56 1.58 1.60 1.62 1.50 1.52 1.54 1.56 

Aobs/*™ 



1.60 1.62 1.64 1.66 



Fig. 4. — Same as Figure|3] but for [O III]-H/3 spectra of the emission line galaxies with S/N>3. The vertical dotted lines in each panel show, from left to 
right, the locations of H,3, [Olll]A4959, [Olll]A5007. 
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TABLE 3 
Star Formation Rates 



ID 



ZHa 



^Ha" 



Lho 



Corrected Lhq 



SFRhc'= 



F24p 



uL, 



u.S fim"^ 



SFRir" 



MODS 11 
MODS 11 
MODS 11 
MODSll 
MODS 12 
MODS 12 
MODS21 
MODS21 
MODS21 
MODS22 
MODS22 
MODS22 
MODS22 
MODS22 
MODS31 
M0DS31 
MODS32 
MODS32 
M0DS41 
MODS41 
MODS42 
MODS42 
MODS42 



■0076 
■0094 
■0274 
■0390 
■0125 
■0255 
■2321 
■2612 
■6131 
■1411 
■2282 
■2658 
■4194 
■5133 
■0033 
■0199 
■0116 
■0153 
■0194 
■0297 
■0112 
■0135 
■0145 



1.6731 ±0.0009 
1.9976 ±0.0002 
2.0808 ±0.0004 
2.3516 ±0.0005 
2.2461 ±0.0003 
2.3979 ±0.0003 
1.7341 ±0.0002 
1.5256 ±0.0011 
2.2678 ± 0.0005 
2.3983 ±0.0002 
2.2087 ± 0.0004 
1.5246 ±0.0035 
2.4862 ± 0.0007 
1.5255±0.0016 
2.4605 ± 0.0003 
1.5681 ±0.0020 
1.4878 ±0.0005 
2. 1865 ±0.0003 
2.2425 ± 0.0007 
1.6855 ±0.0003 
2.1231 ±0.0002 
2.3049 ± 0.0003 
2.0952 ±0.0002 



7.1±2.7 
13.9±1.3 

8.0 ±2.0 
15.1±1.9 

10.0 ±0.9 
8.7±1.1 

8.1 ±0.9 
7.7 ±2.2 
1.2±0.3 

15.2±0.9 
7.9±1.5 
17.4±4.8 
8.0 ±2.0 
7.2±4.5 
11.5±1.1 

10.1 ±2.3 
12.4±2.8 
9.6 ±0.9 
11.8±5.3 
8.2±1.6 
16.7±1.3 
9.7±1.0 
7.8±1.5 



1.3±0.5 
4.0 ±0.4 
2.6 ±0.6 

6.5 ±0.8 
3.9 ±0.3 
3.9 ±0.5 
1.7±0.2 
1.1±0.3 
0.5±0.1 

6.8 ±0.4 

2.9 ±0.6 

2.6 ±0.7 
3.9 ±1.0 
1.1±0.7 

5.5 ±0.5 

1.6 ±0.4 

1.7 ±0.4 

3.4 ±0.3 

4.5 ±2.0 
1.6±0.3 

5.6 ±0.4 
4.0 ±0.4 
2.5 ±0.5 



29.4+1" 

°-16.0 

47 5+"-' 

202.0+lJ 
42.8™f 
15.5+f.i 

O+0.8 
^-3.2 
1+33.8 
^-64.2 



31.6+ 



7.1 
100.1+ 



424.0 



-273.9 



4.9+ 



Q+2.2 
-1.2 

48 1+35-' 
to-r_19 3 

43.7+?!:^ 

30 7+^" 

1 9+' ' 
^■^-0.4 
I 1 5+5.4 

6.9+*' 

9 0+^'^ 

17 3+'° 

44 1+9.0 

11 2+"' 

-14.5 



20.O 



144. 1«- 

i55.i;tr 

232 7+57.2 
^-'^■'-109.8 

989.7+'220 
209.9i|f;^ 

7C 0+176.9 

'^■°-38.8 

38.2+^5% 

490.6;^«f 
3.8+|0 ■ 
2077 4+408.2 

23.9+|0/ 
235.6!^^'/ 
214.2+lf,o 
^oq 1+251.9 

150.6!^|j 
56 3+2«-3 
43 g+ltv 



60.3 



-21.1 
1+43.9 



54.6!5« 
9SAfd 



130.1 ±5.5 


2.5±0.1 


114.6±4.9 


106.8 ±3.9 


3.3±0.1 


165.6±6.1 


106.4 ±3.9 


4.0 ±0.1 


219.8±8.1 


127.2 ±3.9 


6.4 ±0.2 


599.4 ±18.5 


282.1 ±3.9 


12.7 ±0.2 


1581.8 ±22.0 


9.7±9.1 


0.5 ±0.5 


10.2±9.6 


26.6±3.9 


0.6±0.1 


14.6 ±2.2 


99.0 ±3.9 


1.7±0.1 


77.7±3.1 


12.6±3.9 


0.6 ±0.2 


11.5±3.6 


148.3 ±3.9 


7.2 ±0.2 


822.1 ±21.7 


12.7±3.9 


0.5 ±0.2 


10.2±3.2 


I81.1±3.9 


2.9±0.1 


165.6±3.6 


13.4±3.9 


0.7 ±0.2 


20.1 ±5.9 


48.7 ±3.9 


0.8±0.1 


33.1 ±2.7 


113.3±3.9 


6.1 ±0.2 


717.0 ±24.8 


13.8±3.9 


0.2±0.1 


6.6±1.9 


54.0±3.9 


0.8±0.1 


35.9 ±2.6 


71.6±3.9 


2.8 ±0.2 


129.9±7.1 


118.9±3.9 


4.9 ±0.2 


378.8 ±12.5 


123.6±3.9 


2.5±0.1 


114.6±3.6 


19.8±7.3 


0.7 ±0.3 


17.4 ±6.4 


78.1±3.9 


3.7±0.2 


236.5 ±11.9 


32.0±5.2 


1.2±0.2 


33.1 ±5.4 



" Spectroscopic redshift measured from the Hq emission line 

'' Observed flux of the Ha emission line after correcting slit loss, in units of 10"" ergs s^'cm"- 

"^ Observed luminosity of the Ha emission line, in units of 10^^ ergs s"' 

■1 Extinction coirected luminosity of the Ha emission line, in units of lO"*' ergs s"' 

■^ SFR derived from extinction corrected luminosity of Ha emission using Equation|4] in units of MQyr" 

' Observed flux density of MIPS 24 fiva, in units of ^Jy. The fluxes of the objects undetected in the public catalog 

image. 

° Rest-8/.tm luminosity inferred from 24 fiva flux with /:-correction bv lCharv & Elbazl 120011) template, in units of 

•^ SFR derived from MIPS 24 fim flux, in units of M0yr"' . 



are measured from the public 
10*'' ergs s"'. 



expected redshift range of sBzK galaxies, it is the most plausi- 
ble candidate for a strong emission line. The identification of 
Ha is confirmed in all cases when other emission lines, such 
as [Nii]AA6548,6583, [Oiii]AA4959,5007, and H/3 are also 
detected at the expected wavelength. Otherwise, we check the 
consistency of the redshift of the Ha emission line with that 
of the SED model fitting using broad-band data (see 33.4b . 

Although one of the conventional methods to estimate ran- 
dom error is to obtain the rms of the scattering of background 
counts at both sides of the object, the slit length for each ob- 
ject is not long enough. The readout noise is suppressed down 
to 16 e~ by non-destructive readout mode, so that the read- 
out noise is about 1/5 of the background noise. Because the 
readout noise and the photon noise of objects are negligible, 
we estimated the background noise by assuming that the ran- 
dom error is dominated by sky-background noise, which fol- 
lows Poisson statistics. In parallel to the reduction of the ob- 
ject spectra, sky -background spectra were reduced without the 
background subtraction, and the Poisson error of each pixel 
is estimated from them. We found that the Poisson error is 
systematically lager than the rms error by a factor of ^ 1.5, 
though the scatter is large. The rms errors can be underesti- 
mated, because counts in the two-dimensional spectra corre- 
late with the neighboring pixels due to the sub-pixel shift in 
distortion correction and wavelength calibration. 

To measure the flux and ce ntral wavelen gth of the emis- 
sion lines, we used SPECFIT (lKrissl [l99?). which fits mul- 
tiple emission lines and continuum simultaneously. Gaus- 



sian profiles were used for emission lines. The continuum 
was assumed to be a linear function of wavelength. The flux 
of [Nll]A6548 was fixed to be 1/3 of [Nll]A6583 flux. If 
H/3 and/or [Olll]AA4959,5007 emission Hnes were detected 
on the two-dimensional spectra, we measured the flux in the 
same way. The random errors of the fluxes were estimated 
from the best-fit parameters of the model fitting based on the 
sky-background error described above. The emission lines 
detected with S/N>3 are listed in Table [T] and the spectra are 
shown in Figures[3]and|4] The measured flux and the redshift 
derived from the central wavelength of Ha emission line are 
shown in Table [3] 

It is difficult to accurately evaluate systematic flux error 
since part of the flux from the object will be lost by a slit mask; 
the systematic error due to the flux loss is not negligible if we 
measure the total flux of emission lines from slit spectroscopy. 
The flux loss can be caused by a variety of things, for exam- 
ple, a larger size of the object than the slit width, imperfect 
alignment of the slit mask to the target, and misalignment of 
the dithering direction along the slit. 

First we estimated the flux loss for a point source. When 
a Gaussian profile with FWHM = 0."6, which is the typical 
seeing during the observations, and a slit width = 0."8 are as- 
sumed, fluxes passing through the slit are 88.4% of the in- 
tegrated flux; the aperture correction factor is ^ 1.13. We 
applied this factor to the point sources observed. We found 
that, however, the fluxes of the point sources were still under- 
estimated by up to ^ 20 %, even if we also correct the effect 
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of the seeing size variation and pointing error, which are mea- 
sured from alignment stars. This fact indicates that imperfect 
alignment of dithering and/or short-time variation of the see- 
ing size during the observations can also add to the systematic 
error 

On the other hand, it is expected that the flux loss for ex- 
tended sources varies with the size of the sources. We cal- 
culated aperture correction factors from the ratio of total flux 
to that expected to come in through slits, assuming that the 
regions emitting stellar continuum and Ha had the same spa- 
tial distribution. The flux through a slit was measured on the 
^5-band image with a rectangle of the same aperture size as 
for extracting one-dimensional spectrum and the same posi- 
tion angle of the slit, while the total flux was obtained from 
MAG_AUTO of SEXTRACTOR. The alignment error of the 
slit was estimated from the alignment stars on the slit mask. 
We found the median of the aperture correction factors to be 
^1.9 with little dependence on the size of the objects. The 
rms scatter of the factor among the objects is ±0.5. We ap- 
ply the median value to all objects to correct the flux loss. 
We checked the above aperture correction method by com- 
paring the continuum flux derived from the corrected spec- 
trum with total magnitude, which was measured from Ks or 
//-band imaging data according to the redshift of each galaxy. 
We used 7 galaxies in the sample whose stellar continuum 
around the Ha emission line was detected with sufficient S/N. 
We applied the aperture correction to the continuum flux, then 
contrasted the corrected flux with the total magnitude. The 
comparison shows that both results are consistent and the un- 
certainty in the aperture correction is within ~ 20%. 

We need to note that the assumptions on the light distribu- 
tion of Ha emission line would also cause systematic errors. 
The Ha and continuum emitting regions would not always be 
similar for star-forming galaxies at z ~ 2. Studies with the 
integral-field spectrograph show that Ha emission d oes not 
trace the stellar continuum in some galaxies ( Swinba nk et alJ 
I2OOI l2006t iKriek et allllOOTl) . However, this kind of check 
on our slit-spectroscopy sample is difficult, because both the 
continuum and emission lines are so faint that the light distri- 
bution along the slit is not able to be measured. As a lower 
limit of the aperture correction, when the Ha emitting region 
is point like, the flux becomes smaller by ~ 40%. 

3.3. Stacking Analysis 

The flux ratios of rest-optical emission lines such as 
[Nll]A6583/Ha, [Olll]A5007/H/?, and Ha/H/3 provide 
us information on physical properties of ionized regions 
such as the energy sources for ionization ([NllVH a and 
rOlll 1/ H/3; [Bal dwin et al. 198 1 ; Veilleux & Ost erbrockl 
[l987t iKewley et al.1 120011: .KauffmanneLaLl l2003h . die 
amount of dust ext inctio n (Balmer decrement, e.g., 
1989), 



Ha/H/3; lOster brock 1989), and gas-phase metalhcity 
([N Il 1 /Ha: IStorchi -Bergmann et al.l 1 19941: iRaimann et all 
120001: iDenicoloetal.; 2002; Pettini & Pa^|2004D. However, 
except for Ha, the emission lines of the present sample are so 
faint that we can not evaluate the line ratios of an individual 
galaxy. Therefore, we stacked the calibrated spectra in the 
[Nil] -Ha, and [Olll]-H/3 wavelength ranges to obtain high 
S/N spectra, and discuss the averaged physical properties. 
The stacked spectrum F\^ was calculated with 



„ Fnaji 



(1) 



where F\^j,{Xe) and Fhum are the spectrum and the total Ha 
flux of nth object, respectively. The spectra were re-sampled 
in the rest-frame wavelength scale (Xg). The fluxes were 
scaled with total Ha flux, though the variety of fluxes among 
the samples are not large. Because our main purpose of stack- 
ing is to examine the average physical properties of galax- 
ies from their line ratio, weighted mean (by either inverse 
square of noise or luminosity) is not applied; the line ratios 
of weighted mean flux could be biased toward strong emis- 
sion line galaxies. 

We selected the objects whose spectra covered the wave- 
lengths of redshifted Ha, [Nil], H/3 and [Olll]. The emis- 
sion lines of the objects which had at least H/3 or [O III] with 
S/N > 3 were stacked for the measurements of the line ra- 
tios of [N Il]/Ha and [O lll]/H/3. As a result, 9 galaxies were 
used for the stacking. The stacked Ha- [Nil] and H/3-[Olll] 
spectra are shown in Figure |5] The emission lines of the 
stacked spectra were fitted with a combination of the same 
Ga ussia n profiles and linear continuum measurements used 
in 33.21 The best-fit spectra are also shown in Figure |5] We 
also made stacked spectra including the objects with neither 
H/3 nor [O III] detection, but the resultant line ratio was not 
changed. 

To evaluate the effect of galaxies which would happen to 
have a significantly different line ratio from the others, we 
used a bootstrap method. We built 1000 bootstrap sets, whose 
galaxies were randomly selected, allowing duplicate sampling 
from all galaxies used in the stacking, and stacked the spec- 
tra of the galaxies in each set to derive the line ratios. We 
regard the 16% and 84% quantile of the distribution as Icr 
uncertainty in line ratio. 

3.4. SED fitting 

To evaluate the properties of stellar population such as stel- 
lar mass, amount of dust extinction (Ay), and age, we use 
multiband photometric data from the MODS JHKs images 
and publicly available images. The pu blicly availab le im- 
ages we use are: {/(KPNO/MOSAIC; ICapak et alj ^004), 
BViz (ACS/HST; GO ODS Treasury Program Version 1.0; 
iGiavaUsco et"ani2004l) . 3.6 /xm, and 4.5 ^im (IRAC/SST; DR2 
of the Spitzer Legacy Science program for GOODS; Dick- 
inson et al. in preparation). These images are deep enough 
to detect our galaxies with sufficient S/N and to obtain the 
properties of stellar population. The flux of the galaxies were 
obtained with S/N = 8 for > 80% of the objects in all bands, 
except for B (60 %) and U (24 %) bands. The aperture size for 
photometry is set to 1 ."5, which is approximately 2 x FWHM 
of the PSF size in K^ band. Aperture corrections were applied 
to each band properly (see .Kaiisawa et al...2009.) . The con- 
tribution of Ha emission to the flux in K^ or H band, which 
depends on redshift, was subtracted. The corrections are 0.01 
- 0.34 mag with a median of ^ 0.08 mag. 

We fitted the photometry data with stellar popula- 
tion SED models. The SED libraries of GALAXEV 
jBruzual & Chariot 2003) are used for the fitting. We adopted 
the librar i es wit h metallicities of I.OZ0, 0.42©, 0.2Zq and a 
IChabried (120031) IMF. Exponentially declining star formation 
histories (r-model) are assumed: SFR ex e~'^^ with r = 0.01, 
0.02, 0.05, 0.10, 0.20, .50, 1.0, 2. 0, and 5.0 Gyr. The star- 
burst attenuation law of ICalzetti et a l. (2000) was applied to 
the model with 0.1 steps from Ay = 0.0 to 3.0. 

Estimations of the stellar m ass are affected by the adopted 
IMF. Many authors employ a 'Salpeted (Il955h IMF, which is 
known to overestimate stellar mass due to the steep faint-end 
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Fig . 5 . — Same as Figure[3] but for the stacked [O III]-H/3 (left) and the stacked [N n]-Ho (right) spectra. The objects which have [O III] and/or H/3 emission 
are stacked for both spectra. 



slope of the IMF. In addition to the'Chabrier (2003) IMF, we 
fitted the SEDs to the model with the .Salpeter. (.1955 ) IMF 
and measured stellar ma ss for compariso n. For our sample, 
the stellar mass with the Salpeter" ("1955") IMF is larger by a 
factor of - 1 .7 than that with the Chabrier (2003) IMF 

Recent evolutionary population synthesis models, for ex- 
ample Ntoaston (2005), which take account of contributions 
of thermal-pulsating asymptotic giant branch stars to SEDs, 
result in changes i n the in ferred SED properties. We also fit- 
ted the SEDs with lMarast on (2005) model and found that the 
model made the stellar masses smaller by a factor of ^ 1 .6. 

HYPERZ (Bolzonellaet al. 2000) was used to determine 
best-fit SED models. The redshift for the fitting was fixed 
to that derived from the Ha emission line. We derived to- 
tal stellar masses of galaxies, scaling the stellar masses and 
SEDs with the total magnitude (MAG_AUTO) of the K^ band. 
The parameter matrix, which gives a y^ value for each grid 
of SED models, was used to evaluate the uncertainty in the 
model parameters. The uncertainty is obtained from the prob- 
ability distributions of the parameters. For example, the prob- 
ability distribution as a function of stellar mass (f (M*)) was 
calculated from a minimum x^ value for each stellar mass 
ix^iMi,)). We determined a 68% confidence level of the dis- 
tribution as Ict error of the parameter. The stellar population 
parameters and the Ict errors derived from the SED fitting are 
listed in Table |4] Some examples of the best fit SEDs and 
the probability distributions of stellar mass. Ay, and age are 
shown in Figure |6] 

4. RESULTS 

4. 1 . Properties of Ha Emission Lines 

Ha emission lines are detected from 23 out of the 37 galax- 
ies (62 %). The detection rate does not depend on the mid-IR 



flux. The rate for sBzK-MIPS and sBzK-non-MIPS galaxies 
are 11/18 (61 %) and 12/19 (63 %), respectively. 

The spectroscopic redshifts are consis tent with theu" photo - 
metric redshifts in the MODS catalog ( Kajisawa et al.ll2009h : 
the median value of the relative errors is {dz) ~ -0.026 with 
a scatter of 0.035, where dz = (Zphot-ZHa)/(l +ZHa)- The 
spectroscopic redshifts are newly determined for 10 galax- 
ies. Out of the remainin g 13 galaxies with spectroscopic 



redshifts in the literature ( CohenI 120011; ICohen et all 



)henl 

Dawson etal. 2001; Wir th et'alThoO^ iCowie et al. 



i Treuet all 12 005; Chapma n et ail 120051; iReddvet al.1 



200C; 



2004 



2006bt 



iBarger et al . 2008), two galaxies have slightly different spec- 
troscopic redshift. MODS32-01 16 has zna = 1 -488 (ziiteratm-e = 
1.315) and MODS2 1-2321 has zua = 1-734 (zuterature = 1-879). 

The left panel of Figure [T] shows the apparent Ha lumi- 
nosity and the redshift distribution of the Ha-detected sBzK 
galaxies. The galaxies lie at redshifts between 1.49 and 2.49 
with a median of 2.12. The redshift distribution is con- 
sistent with that expected for BzK galaxies (1-4 < z < 2.5; 
iDaddi et al.l2004l) . The lack of galaxies at 1 .8 < z < 2.0 is due 
to an atmospheric absorption band between H and K bands. 

The luminosity distances of these objects are calculated 
with their redshifts measured using the Ha emission lines. 
The fluxes of Ha emission lines (Fna) range between 1 and 
17 xlO~'^ergss"' cm"^, which corresponds to an Ha lumi- 
nosity (Lho) of 0.5 to 7.0 X lO'^^ergss"'. Fhu and Lho are 
listed in Table [3] 

Possible reasons for non-detections of emission lines from 
the remaining 14 galaxies are; (1) their Ha emission lines are 
too faint for the present observations; (2) the redshift of those 



galaxies are out of the 1 .2 



<7< 



2.5 range, where the grism 



(HK500) we used covers the Ha wavelength; (3) the galax- 
ies lie at the redshift where Ha emission line falls on strong 
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Fig. 6. — Examples of best fit SEDs (top) and probability distribution against stellar mass (middle) and age versus the amount of extinction (bottom). Left 
and Right panels show MODS12-1025 (non-MIPS) and MODS22-2658 (MIPS), respectively. In the middle diagrams, the dashed lines show the range of 68% 
confidence intervals (Icr error), and the thick dashed line shows the best-fit stellar mass. 



sky emission lines or atmospheric absorption bands. Because 
the photometric redshifts and the errors of the galaxies with 
emission lines undetected are at 1 .2 < z < 2.5, the second pos- 
sibility is unlikely except for MODS2 1-2047, whose znterature 
is 0.999. In order to check the third possibility, we calculate 
the detection limit of Ha emission as a function of redshift 
using the sky background noise and the response spectrum. A 
Gaussian profile is assumed for Ha with FWHM = 540km s"', 
which corresponds to the wavelength resolution of the HK500 
grism with a 0."8 slit width. Provided that Ha emitters have 
a uniform redshift distribution at 1.4 < z < 2.5 and an appar- 
ent Ha luminosity of 1.8 x 10^^ ergs s~', which corresponds 
to the 84 % quantiles of the distribution of apparent Ha lu- 



minosity, the expected detection rate of the Ha emission line 
with S/N > 4 is 87.7 %, which corresponds to ^ 32 objects: 
objects at 1-8 < z < 2.0 are rarely detected. Therefore, for 
5 of the observed samples. Ha emission lines are undetected 
due to the strong sky emission lines or atmospheric absorp- 
tion bands, while the remaining 9 galaxies have smaller SFR 
and/or larger attenuation than the galaxies with emission line 
detection. 

4.2. [O m]X5007/Hl3 and [Nu]\6583/Ha Line Ratios 

Diagnostic diagrams based on optical emission-line ratios 
are commonly used to determine the energy source of ioniza- 
tion. Possible sources are massive stars, AGNs and energetic 
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Fig. 7. — Luminosity of Ha emission lines and redshift distribution of the Ha-detected sBzK galaxies before (left) and after (right) extinction correction. Red 
and black symbols show sBzK-MIPS and sBzK-non-MlPS galaxies, respectively. In the right figure, small symbols show the results of the equal-extinction case 
correction to the stellar continuum. 



shocks (iBaldwin et al.ll98Ulwileux & Osterbrockll987h . In 
this paper, we use Ha luminosity as an indicator of SFR, as- 
suming that its energy source is massive stars. We excluded 
AGN candidates from the sample using multi-wavelength data 
(see ^2.1b . In addition to these diagnostics, we use optical 
emission-line ratios to determine the energy sources of the 
present emission line objects. 

In Figure [H we plotted [Olll]A5007/H/3 and 
[Nii]A6583/Hq; line ratios of the stacked spectra and 
the galaxy with the four emission lines detected (MODS31- 



0033). Results for ~ 347000 galaxies from SDSS DR6 
jAdelm an-McCarthv et alj 120081) and local emission line 
galaxies ftMoustakas & KennicuttI l2006l) are also shown. A 
classific ation hne empirically d etermined from the SDSS 
objects (Kauffmannet al. 2003) and a theoretically deter- 
mined classification line (Kewlev et al. 2001) are also shown. 
Emission-line galaxies located below and to the left of the 
lines are expected to have massive stars for their energy 
sources. The distribution of local SDSS galaxies exhibits the 
two expected sequences: the populations below (star-forming 
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Fig. 8. — [Nn]-Ha and [OlII]-H/3 line ratios. The black diamond shows 
the line ratios of the stacked spectra, while the white diamond is the stacked 
spectra with MODS41-0194 (a possible AGN candidate; see text) excluded. 
The open square is MODS31-0033, in which the four emission lines are de- 
tected. The contours in gray scale show the distribution of ~ 347000 objects 
from SDSS DR6, while the black dots show local starburst galaxies analyzed 
by Moustakas & Kennicutt (2006). The solid line is the empirical ly deter- 
mined classification line derived for the SDSS objects by Kauffm ann et al.l 
Q003). Objects below and to the left of this line are normal starburst galax- 
ies. The dashed line is the classification boundary theoretically determined 
bv lKewlev et alj 120011) . 

galaxies) and above (AGNs) the solid line. 

The emission line of the stacked sample place it on the the- 
oretical line, while MODS3 1-0033 is located to the lower left 
of the empirical line. Our results are located above the se- 
quence of local star-forming galaxies, closer to the AGN se- 
quence. In fact, a similar offset is seen in star-forming galax- 
ies at z ~ 2 d&b et al. 2006a; Kiiek et al. 2007). One pos- 
sible reason for the offset is a higher electron density and a 
harder ionizing s pectrum than for local star-forming galaxies 
(lErb et alJl2006al) . 

We note that AGNs can be also identified based on [O III]- 
Ha line ratios for objects without an H/3 detection. AGNs 
typically have [O lll]/H/3 > 3 on Figure |8] On the other 
hand since the Ha/H/3 line ratio must be ^ 3 (the B aimer 
decrement) or larger (if there is a significant dust extinction) 
the [O III] -Ha line ratio of an AGN candidate will be larger 
than 1. Only one object (MODS41-0194) has [Olll]-Ha > 1 
(1.67). When we exclude MODS41-0194 from the stacking, 
the stacked line ratios slightly move below the empirical line 
(Figure [8]|. We regard this object as a possible candidate of 
AGN in the later discussion. 

4.3. Extinction 

In order to evaluate SFR from intrinsic Ha luminosity, the 
extinction due to interstellar dust needs to be corrected. More- 
over, the samples are expected to have large amounts of dust, 
because about half of the present galaxies are luminous in- 
frared galaxies, in which strong infrared emission is thought 
to originate from dust thermal emission. 

The extinction of emission lines in star- forming regions can 
be estimated with the reddening of the Balmer decrement 
through an extinction curve. The intrinsic line ratios of the 
Balmer emission lines are nearly constant with respect to elec- 
tron temperature of ionized gas. The amount of reddening is 



Fig. 9. — Comparison of dust extinction derived from an SED fit to that 
derived from the Balmer decrement. The black diamond shows the stacked 
spectrum. Red and black squares show sBzK-MIPS and sBzK-non-MIPS 
galaxies, respectively. Squares with error bars are the galaxies with an H/3 
emission line detected (MODSll-0390, MODS12-0255, MODS22-1411, 
MODS31-0033), while those with upper an'ows are the galaxies with an H/3 
upper limit (3a-). The hatched region shows the 3(t detection limit of H/9 
estimated from Ho flux. The height of the region is calculated from the 
68% distribution of the Ha fluxes. The dotted and dashed lines illustrate the 
equal-extinction case and the relation given bv .Calzett i (200 IJ (Equation[3}, 
respectively. 

estimated from the ratio of the attenuated-to-intrinsic Balmer 
fine ratio, 



E(B-V)„ 



2.5 



7^Hc 



log- 



/F, 



^Hfi 



KXho) - KXhp) Fo^"" /Fo 



\h0 ' 



(2) 



where "i" and "o" denote intrinsic and observed flux, re- 
spectively, an d k{X} is an extinction curve. We assume 
ICalzettietall ([2000) law for k(\). 

The Balmer decrement method can only be applied to the 
14 galaxies at z > 1.9 in the present sample, because the H/3 
wavelength is not covered by the HK500 grism for galax- 
ies at lower redshift. Both Ha and H/3 are detected in four 
galaxies (MODSll-0390, MODS12-0255, MODS22-1411, 
MODS3 1-0033), so that E(B-V),.^, is estimated using the 
above equation. For the rest of the sample, lower limits on 
the E(B-V)g3s are evaluated using H/3 emission-line upper 
limits. The resultant £'(B-y)gas values are > 0.5. We also 
used the stacked spectrum to measure the average amount of 
reddening and found it to be E(B -V)gas = ^■'^^^)02- 

We estimate the amount of reddening for the ionized gas 
from the SED fitting and apply these extinction corrections to 
the Ha luminosities for all of the galaxies. Figure |9] shows 
the comparison between the amounts of reddening derived 
from the Balmer decrement and SED fitting. Our result shows 
that the gas is more attenuated than the stars. From SED fit- 
ting, the amount of reddening for stars in the samples used 
for stacking are found to be E(B -V)ste\iir = 0.20^oQg, thus 



, 26+''-42 

^■^"-0.10- 



It should be noted that 



£(B-V),stelto/£(B-V)gas 

the galaxies which we used to determine the relation suffer 
relatively less attenuation since we would not be able to de- 
tect H/3 emission in the more highly attenuated galaxies. The 
hatched region in the figure also shows the lower limit of 
£'(B-y)gas corresponding to the 3a detection limit of H/3. 
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The height of the region calculated from the 68% distribution 
of the observed Ha fluxes. 

From studies of local starburst galaxi es, it is k nown that 
ionized gas is more attenuated than stars (ICalzetti et al. 2000; 
lCalzettill2001h . The difference is thought to be caused by 
changes in the ge ometry of the ionized gas and stars in a typi- 
cal galaxy. ICalze tti (2001) give the relation between the color 
excess observed for gas and stars: 



£(B-V)steiiai- = 0.44£(B-y)g 



(3) 



This relation is shown with the dashed line in Figure |9] Our 
stacked result is consistent with the equation within errors. On 
the other hand, H/3 detected galaxies are also consistent with 
the case that extinction for gas and stars is equal (defined as 
"equal extinction", hereafter), though these galaxies would lie 
at the lower end of reddening as noted above. We, therefore, 
employ Equation[3]as an averaged nature of dust extinction to 
correct for the effect of dust extinction on the emission lines 
and compare these with that calculated assuming equal ex- 
tinction. We discuss the reliability of th e application of these 
extinction correction further in the ^5.11 

The right panel of Figure|7]shows the luminosity of the Ha 
emission lines of the galaxies after the extinction is corrected. 
The extinction correction factor of our samples ranges from 
1 .2 to 87, resulting in a wider spread in extinction corrected 
Ha luminosities (^ 0.6 dex for the CalzettI (120011) relation 
and ^ 0.4 dex for equal extinction) than in the apparent Ha 
luminosities (~ 0.3 dex). This result indicates that proper ex- 
tinction correction is important to reliably calculate the intrin- 
sic Ha luminosity. 

4.4. Star Formation Rate and Stellar Mass 

To convert the extinction-corrected Ha luminosity to SFR, 
we derive a relation between them using stella r popula- 
tion synthesis models, following 'KennicuttI £ 9981). We em- 
ploy GALAXE V models (Bruzual & Charlotl l2003h with the 
IChabried (l2003h IMF and the same parameters as used in ^ 
to calculate ionizing photons (A^Ly)- N^y i s converted to L(H a) 
based on a case-B recombination model ftOsterbrockl 19891) . If 
a constant SFR model and solar metallicity are assumed, the 
L{Ha)-SFR relation with an age older than 20 Myr is given 

by, 

SFR(Ha)(MQ yr"') = 4.9 x lO-'^^L(Ha)(eTg s"'). (4) 

Using this equation, we derived the SFR with extinction- 
corrected Ha luminosity. The result is listed in Table |3l The 
SFR estimate is larger by a factor of ^ 1.5, if the iSalpeterl 
([1955) IMF is adopted. 

The stellar masses derived from the SED fitting are com- 
pared with the SFRs calculated using both extinction cor- 
rection appro aches in Figure [TOl For the extinction correc- 
tion using the ICalzettil (120011) relation, the SFRs range from 
10 to several hundreds Moyr"'. Some massive galaxies 
(~ 1O"M0) and a significant fraction of low mass galaxies 
show SFR larger than lOOMoyr"'. As a result, there is no 
clear correlation between SFR and stellar mass. The present 
sample contains many sBzK-MIPS galaxies, whose SFRs are 
expected to be high. Although sBzK-MIPS galaxies lie at the 
massive and high SFR end of the distribution, there still ex- 
ist low-mass and high-SFR galaxies among sBzK-non-MIPS 
galaxies. 

The specific SFR (SSFR), which is the SFR divided by stel- 
lar mass, is one of the indicators showing the level of star for- 



mation; higher SSFR indicates that a galaxy is accumulating 
large stellar mass compared to its own stellar mass on a short 
timescale. In Figure [TO] SSFR is depicted by dashed lines. 
SSFR does not depend on MIPS detection but on stellar mass. 
There exist actively star-forming galaxies (SSFR > lO'^^yr"') 
among the less-massive galaxies, while the SSFR of massive 
galaxies (> lO'^M©) is nearly constant (-- 2 x 10"'^ yr"'). It 
can imply that some low mass galaxies at z ~ 2 are actively 
forming stars. 

On the other hand, for the SFRs corrected assuming the 
equal-extinction case, the SFRs and SSFRs of the high- 
SSFR galaxies become lower, but they still have higher SSFR 
when compared to other galaxies. The SSFR of the original 
high-SSFR galaxies become SSFR > 3.5 x 10"Vr"'. Conse- 
quently, there is still no correlation between them. The distri- 
bution of SFRs is nearly constant {^ 3OM0 yr"'). 

The possible AGN candidate (MODS41-0194; denoted by 
a open square) shows lower SFR than other galaxies with sim- 
ilar stellar masses. We note that this object does not change 
our findings of no correlation between SFR and stellar mass 
and low mass galaxies with active star formation. 

The inverse of SSFR (SSFR age) gives the estimated age 
for a galaxy to build up all of its stellar mass with the current 
SFR. On the contrary, mass-weighted age, which is calculated 
by integrating the ages of stars weighted by SFR at that time, 
is an indicator of the average age of the stars in a galaxy. The 
mass-weighted age is given from integration of the star for- 
mation history of galaxy as 



['^(t')x{t-t')dt' 

Tit) = -^Q ; — - — 



(5) 



where ^(f) is the SFR at age t. Comparing the two ages, 
we can roughly check consistency between the SED fitting 
result and the SFR estimated from extinction corrected Ha 
luminosity. The comparison between the SSFR age for the 
both cases of extinction correction and mass-weighted age of 
our sample is shown in Figure [TT] For extinction correction 
with the ICalzettil (1200 ll) relation, the SSFR ages agree with 
the mass-weighted ages to within 1 a error; a galaxy with a 
higher SSFR has a younger stellar population. This fact sug- 
gests that the majority of the stellar mass for the galaxies in 
our emission line samples were formed in a recent starburst. 
On the other hand, for the extinction correction assuming the 
equal extinction case, the ages of stars are young compared to 
the SFR especially for high-SSFR galaxies. In this case the 
galaxies must have experienced much more active star forma- 
tion in the recent past and terminated the star formation in a 
very short timescale (< 10 Myr). 

5. DISCUSSION 

5.1. Extinction Correction and Star Formation Rate 

The median value of the reddening derived from the SED 
fitting is £'(B-y).steiiar = 0.42, 0.21 (Ay = 1.7, 0.85) for sBzK- 
MIPS and sBzK-non-MIPS galaxies, respectively. Both 
values are relatively larger than those for local galaxies 
(e.g., (Ay)^ ,^;.,„ ^ 0.43 for r-band selected SDSS galaxies; 
iKauffmann et al. 2003). The extinction corrections are key to 
the discussion on the luminosities of Ha emission line, and 
thus the SFRs and SSFRs; larger extinction correction makes 
extinction-corrected SFR and SSFR larger In this section, 
we discuss the reliability of our extinction correction through 
comparisons of the Ha SFRs with other SFR indicators, and 



14 



Yoshikawa et al. 



Calzetti relation 



LO"'i/r" 



4^ 



LO-'i/r-', 




10"'°i/r-' ,'' 



10'" 



1.5 « 
C/3 



0.5 
0.0 



10*- 



equal extinction 



10"' yr" 



L0"»i/r" 



LO"°i/r- 



P-T^ 







*>:*. 



10"'°i/r-' ,-' 



10'" 



1.5 « 
Co 



Fig. 10. — Ha luminosity with extinction collected using the Calzetti relation (/e/f) and the equal-extinction case (right) as a function of stellar mass derived 
from SED fit ting Red and black squares show sBzK-MIPS and sBzK-non-MIPS galaxies, respectively. MODS41-0194, which is the possible AGN candidate 
discussed in JI4.2I is denoted by an open square. SFR//a (the right axis) is calculated by Equation. |4] The dashed lines indicate SSFRs. 
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Fig . 11 . — SSFR age (SSFR ' ) with extinction corrected using the Calz etti relation (left) and the equal-extinction case (right) as a function of mass weighted 
age. The dashed line shows the equal age. Symbols are same as Figure fTol 



the reason why the higher SSFR galaxies are found by our 
observation. 

5.1.1. UV Luminosity 

The UV continuum is used as a SFR indicator, because 
the stellar continuum at 1500-2800 A is sensitive to emission 
from massive stars. We compare the SFRs with those esti- 
mated from the UV continuum. We calculate a relation be- 
tween UV continuum at A = 1500 A, L,y.i5oo, and SFR as for 
Equation]?] Li/.isoo is measured from the S ED te mplates of 
GAL AXE V with the same parameters as in 34.41 with the B- 
band filter of ACS/HST redshifted to z = 2. The relation is 
given by, 

SFR(UV)(M0 yr-') = 7.50 x 10-2X,i5oo(erg s"' Hz"'). (6) 



For our samples, Liy.1500 is measured with (pl'.'S in B-band 
with an aperture correction obtained from the /Tj-band total 
magnitude. The dust extinction is corrected with that derived 
from the SED fit. 

The comparison between SFR(Hq;) and SFR(UV) is shown 
in Figure 112] For the equal-extinction case, the SFR(Hq;) is 
underestimated for the galaxies with h igher SFR(UV). On the 
other hand, using the ICalzettil (1200 Ih relation the SFR(Ha) 
and SFR(UV) are roughly consistent for a broad range of 
SFR, but are systematically larger by 0.3 dex for both sBzK- 
MIPS and sBzK-non-MlPS galaxies. 

We note that the luminosity ratio of the UV continuum to 
Ha varies with the population of a galaxy. For example, 
because the UV continuum is contaminated by less massive 
stars than stars emitting ionizing photon, galaxies populated 
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Fig. 12. — (left) Comparison of SFRs inferred from Ha luminosity (SFR(Ha)) and UV luminosity (SFR(UV)). Diamonds and squares are high-SSFR 
(SSFR > 10" yr"') and low-SSFR (SSFR < 10 yr" galaxies, respectively. Red symbols are MlPS-selected samples. The SFR(UV) is calculated by Equation 
|6] with UV luminosity, which is measured with ^I'.'S in fl-band with an aperture correction obtained by A^,-band total magnitude. The dust extinction is corrected 
with that derived from SED fit. (right) Same as left, but lCalzettil 120011) relation is not used for extinction connection for Ha luminosity. 



by younger massive stars are more luminous in Ha. For our 
template, the ratio is larger by factor ^ 3 for age< lOMyr 
and becomes constant for age > 100 Myr. In fact, high-SSFR 
galaxies have a mass- weighted age < 100 Myr. If a top-heavy 
IMF is assumed, moreover, the ratio gets larger, so that more 
ionizing photons are emitted than UV continuum. 

5.1.2. MIPS Flux 

In order to examine the reliability of the extinction correc- 
tions further, we compare the SFR(Hq:) estimates with those 
derived from the thermal infrared luminosities (SFR(IR)), be- 
cause the infrared luminosity is an SFR indicator almost un- 
affected by dust extinction. We use the rest-8 fim luminosity 
(i^Li,,fi,^i„) derived from the MIPS 24 /im flux as representa- 
tive of the thermal infrared luminosity. Independently from 
the public catalog, the 24 /im fluxes of the present galaxies 
are measured using the Spitzer MIPS image (Dickinso n et al. 
in pre paration) with the IRAF/DAOPHOT package ( StetsonI 
119921) for both sBzK-MIPS and sBzK-non-MIPS samples. 
Some galaxies have 24/im fluxes severely contaminated by 
nearby bright 24/im sources because of the large PSF of MIPS 
(FWHM ^ 5 "4). In order to reduce the effect of the contam- 
ination, the fluxes are measured by fitting crowded objects 
simultaneously with the PSF measured from isolated point 
sources. The objects are centered using the /riband light dis- 
tribution. The photometric error is estimated on a blank-sky 
region on the image. In order to take into account the resid- 
uals of the fitting, the image from which the fitted models 
are subtracted is measured and the residual is added to the 
error. Then, vL^%^„i is estimated from the spectroscopic red- 
shift of Ha emission line and the observed 24 /im flux densi- 
ties with ^-corrections derived from the luminosity dependent 
SED models of Chary & Elbaz (2001) . The measured MIPS- 
24 /im flux and estimated rest-8 /xm luminosity are also listed 
in Table E] 

The derived I'L^.s.fj.m and Ha luminosities are plotted in 
Figure [T3] The H a luminosities are corrected using both 
the .Calzettil (1200 lb relation (left) and the equal extinction 



case (right) for comparison. The expected relation between 
vLy'i^,,, and Ha luminosity is also shown using a solid line 
in the figure. The relation is derived based on that expected 



between vL, 



v.8/^??? 



and SFR. The vL, 



i/,S/xm 



is converted to to- 



tal infrared luminosity (Lir) in the wavelength range between 
8 and 1000 /xm with the luminosity dependent SED models 
dCharv & Elbaz 2001), and the relation between the total in- 
frared luminosity and SFR, 



SFR(/;?)(Mq yr-') = 3.0 x IO'^^r (ergs"') 



(7) 



10, 



1.2x10-'%r(Lo) 



is used. The relation is derived by divid i ng that in lKennicutJ 
(1998) by 1.5 to convert flie ISalpeted (119551) IMF to the 
ChabrieB (l2003h IMF. As we see in Figure [ill it is sug- 
gested that the extinction correction between the stellar con- 
tinuum and ionized gas depends on the SFR of a galaxy. For 
higher ( l ower) SFR, the Ha luminosity corrected with the 
ICalzettil (120011) relation (equal-extinction case) is consistent 
with vLi,,i^m. We divide the present sample into two groups, 
LIRGs and ULIRGs, at vLy,^^,„, = 2.8 x 10'*'* ergs"', which 
corresponds to Lir = IO'^Lq. For ULIRGs, the systematic 
offset between SFR(IR) and SFR(Ha) (SFR(IR)/SFR(Ha)) 
is 0.31 ( 0.72) dex, with a scatter of 0.47 (0.29) dex for the 
'Calzetti' ('2001') relation (equal-extinction case); although the 
Calzetti (2001) relation results better consistency, SFR(Ha) 
is still underestimated. For LIRGs, the systematic offset is 
-0.42 (0 . 02) de x with a scatter of 0.29 (0.43) dex for the 
ICalzettil (120011) relation (equal-extinction case); using the 
equal-extinction case results in a good agreement between the 
two SFR indicators. 

We need to note that SFR estimated from the rest-8 /im lu- 
minosity has a large uncertainty. Rest-frame luminosity at 
^ 8/im contains a very small fraction of thermal infrared 
emission, so that the Lir is estimated by extrapolating the SED 
to longer wavelengths. The rest-8 /im luminosity could also 
be affected by Polycyclic Aromatic Hydrocarbons (PAHs) 
emission lines, whose energy source is dominated by late 
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Fig. 13. — (/e/f) Comparison of Ha emissions with rest-8 ^m luminosity. Symbols are same as in Figure [12] The 24 /^m fluxes are measured from a Spteer 
MIPS image (Dickinson et al. in preparation) with IRAF/DAOPHOT package ( Stetson 1992) (see text). The solid line shows the expected relation between 
vLyjuj, , and Ha lumin osity (see text). The eiTor bars of vLug/j,,, are estimated from photometric error of 24 /ira flux and the residual of fitting, (right) Same as 
left, but lCalzett] 420011) relation is not used for extinction correction for Ha luminosity. 



short- lived B-type stars rather than O-type stars jPeeters et alj 
|2004). These factors can produce large random error in the 
estimatio n of SFR from the res t-8 /im luminosity. For local 
galaxies. ICharv & Elbazl ( 1200 Ih demonstrated that the la er- 
ror for the LiR estimated from rest-6.7 /im was a factor of two. 
The uncertainty is consistent with the observed scatter around 
the solid line in Figure [T3] 

5.1.3. Radio and Far-Infrared 

A correlation between FIR luminosity and radio continuum 
is known for local galaxies (Condon 1992; Yunet al.]l2001l) . 
Non-thermal radio continuum is mainly emitted by syn- 
chrotron emission from the remnants of supernovae, whose 
frequency is proportional to the number of massive stars, and 
thus the SFR. The relation between these luminosities for lo- 
cal galaxies is given by 

Lir(Lo) = 3.5 X 10-'2Li.4ghz(WHz-') (8) 

with about 40% of uncertainty (lYunet al.ll2001l) . 

We u se a pubhcly availabl e 1.4 GHz radio catalog and 
images dBiggs & IvisonI 120061) . whose la rms noise in the 
GOODS-N region is 5.8 ^Jy beam"' (the beam size = 1.52 x 
1.51 arcsec^). The catalog contains sources with a > 5cr 
detection. Although our Ha-emission-line galaxies are not 
matched with the sources in the catalog, some galaxies are 
marginally detected (< 3a), according to our eye-ball check 
on the 1.4 GHz radio images. Therefore, their 1.4 GHz fluxes 
are lower than 29 fi}y (5 a), and some galaxies possibly have 
~ a few 10 ^Jy. 

If we assume the local radio-lR correlation, we can cal- 
culate an expected Li.4ghz from SFR using Equation [T] and 
[8] We assume a spectral index of (a; fi, oc v'") of the non- 
thermal radio continuum a = 0.8 when we calculate the ex- 
pected 1 .4 GHz flux density. The expected 1 .4 GHz flux esti- 
mated from the MIPS -24 /im flux exceeds the detection limit 
only for one emission line galaxy (MODS 12-0125), even if 
the uncertainty in the estimation of Lir from rest-8 /im lumi- 
nosity is considered. This object is a sBzK-MlPS galaxy and 



has a low SSFR. For this object the SFR(Ha) is also smaller 
than SFR(IR), so that only the rest-8 /im luminosity exceeds 
that at other wavelengths. This infrared excess when compare 
to the radio luminosity can be explained by a warm dust com- 
ponent, which is thought to be cau sed by a compa ct nuclear 
starburst or dust-enshrouded AGN ( Yunet al.l[2001h . 

On the other hand, when we calculate the radio flux from 
the ex tinction-corrected Ha luminosity using the Calzettl 
(12001b relation, the radio emission of 5 emission line 
galaxies (MODS22-1411, MODS21-2612, MODS22-5133, 
MODS 11 -0390, and MODS22-2658) exceed the detection 
limit. The first four galaxies are high-SSFR galaxies. Al- 
though most of these galaxies have marginal detection of 
1.4 GHz flux, the expected radio fluxes ( Li.4ghz '^50- 
150/iJy) are higher than the observed radio flux, even if the 
uncertainty in the radio-lR correlation is taken into account. 
On the other hand, for the Ha luminosity calculated assum- 
ing the equal-extinction case, none of the galaxies exceed 
the detection limit. This would suggest that the amounts 
of exti nctions in the io nized gases are overestimated by us- 
ing the ICalzettil (120011) relation, though the uncertainties in 
their extinctions are large. We also note that, equal extinction 
in turn underestimates the extinction for ULlRGs (see Fig- 
ure [T3]l. These Ha-excess galaxies consist of three ULlRGs 
(MODS 1 1 -0390, MODS22- 141 1 , and MODS22-2658). 

Thermal emission from dust heated by massive stars 
(^ 30K) is observed at sub-millimeter wavelengths. The 
GOODS- N region is a l so ob served with SCUBA/JCMT 

f;0/im; 'Holland elalj [19991) . AzTEC/JCMT (1.1mm; 
Isonetal. 2008), and MAMBO/IRAM (1.2 mm; 
Krevsaetal. 1998) dow n to ^ 2mJy (e.g., Pooe et alj 
2005, 2006), - lmJv(Pereraetal.ir20 08; Chapin et al. 2009), 
and ^ 0.7 mJy dGreve et al.ll2008l) . respectively. One of our 
sample, MODS3 1-0033, is identified with an AzTEC source, 
AzGN27, w hose 1.4 GHz and 2 4 /im counterparts are also 
identified by IChapin et al.l (1200 9^. although this object is not 
detected at both 850 /im and 1 .2 mm wavelengths. 
Interestingly, the 1.1 mm flux (2.3l![ jfjmJy) of MODS31- 
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0033 is consistent with that expected fr om the 24 /im flux 
(1.6mJy), when flie lCharv & Elbazl ( l2001h SED templates are 
assumed, as though the SFR(Hq;) is lower than the SFR(IR). 
None of our samples including MODS3 1-0033 are detected 
at 850 /im, in spite that some of them have SFR(Ha) and/or 
SFR(IR) exceeding several hundred MQyr"'. It is known 
that 850 /im selected galaxies seem to have cooler thermal 
emission than local ULlRGs, which results from the fact that 
the 850 /im flux is lower th an that expected from total IR 
luminosity (Po pe et al] 120 06V Therefore, it is possible that 
850 /im emission is not detected from all ultra-luminous star- 
burst galaxies at z ■^ 2. 

5.1.4. X-ray 

The linear relation between X-ray l uminosity and SFR 
is calibrated with local galaxies (e.g., iGrimm et al.l 120031; 
iRanalli et al.ll2003T) . X-ray emission of star-forming galaxies 
is dominated by high mass X-ray binaries (HMXB), whose 
progenitor masses are > 8Mo and thus they are short-hved. 
The c orrelation of X-ray luminosity with UV (iDaddi et al.l 
120041) and Ha (lErbet alJi2006b) luminosity have also been 
investigated for high-z galaxies. 

Due to the sample selection, the present samples are not 
detected in the Chandra 2 Ms catalog, whose limiting flux 
is 2.5 X 10-1^ erg cm-^s-i in soft X-ray band (0.5-2.0keV), 
which detects the redshifted hard X-ra y fluxes of ob j ects a t 
z ^ 2. According to the calibration bv iGrimm et al.l (l2003h . 
the flux corresponds to SFR(X-ray)= 74— 186Mi7^yr~' for our 
galaxies (the IMF is converted to the iChabrieil (I2003h IMF). 
Even if the errors of the SFRs are taken into account, five 
of our samples (MODS 11 -0390, MODS21-2612, MODS22- 
1411, MODS42-0112, MOD S22-2658) have the extinction 
con-ected SFR(Hq:)s with the ICalzettil (1200 ll) relation larger 
than the upper limit of the SFR(X-ray). The first four galaxies 
are high-SSFR galaxies, and the galaxies except MODS42- 
0112 are also Ha excess galaxies against radio flux. This 
would also suggest th at the am ounts of extinctions are over- 
estimated by usin g the ICalzettil (12001.) relation. We note that 
lErb et al.) (I2006bl) reported the higher Ha luminosities of their 
BM/BX galaxies than the X-ray luminosities. The evolution 
of the first massive stars to HMXB takes '^10- lOOMyr, 
which is comparable to the age of our high-SSFR galaxies. 

5.1.5. Extinction and Age 

Extinction and age could be degenerate in SED fitting, be- 
cause both parameters make the SED redder. For example, 
degeneracy in extinction and age is found in the probability 
distribution for some galaxies; see the bottom panels of Fig- 
ure |6] However, if we limit the age to be larger than 50 Myr 
in the SED fit, for example, then the estimated value of ex- 
tinction is almost the same as that of the original fit, showing 
that the degeneracy has only a marginal effect on the estimate 
of extinction. The choice of star-formation history applied to 
the SED model also affects the estimate. We use r-models 
in the fitting, but the constant star-formation (CSF) model 
tends to result in larger attenuation. In fact, if we apply the 
CSF model to our galaxies, the E(B-V)s for about half of the 
galaxies with age > 1 Gyr become larger by AE(B-V) > 0.1 
than that with the r-model, while the E(B-V)s are unchanged 
for younger galaxies. In short, the CSF model gives similar 
or larger best-fit extinction than the r-model in general. Al- 
though a detailed analysis of the rest-optical continuum would 
reduce further the possible range of solutions for best-fit pa- 



rameters dKriek et al.ll2006l) . it is difficult to measure the con- 
tinuum from the present spectrum data. 

5. 1 .6. Comparison with Other Studies 

Our result sh ows no corr e lation b etween stellar mass and 
SFR. However, Dad di et all (l2007ai) demonstrated that BzK- 
selected star-forming galaxies have a tight correlation within 
0.2 dex using SFR derived from UV and mid-I R. We com- 
pare ou r result for SFR and stellar mass with the 'Dad di et alj 
( I2007al) relation (Figure [Hi. For the ICalze tti (2001) rela- 
tion, significant fractions of sBzK galaxies (4/18 and 4/19 for 
sBzK-MIPS and sBzK-non-MIPS, respectively) deviate from 
the correlation. Massive galaxies (> IO^'Mq) in our sam- 
ple are distributed around the relation given by iDaddi et aP 
(l2007al) . while some of the less massive galaxies (< lO'^M©) 
have much higher SFR than that expected from the correla- 
tion. Even for the equal-extinction case, which makes the 
SFRs of the higher SFR galaxies lower, there is sti ll no corre- 
lation b etween stellar mass and SFR like that of IDaddi et alj 
(lIOOTal) . 

One of the reasons for the discrepancy could be the dif- 
ferent extinction estimates used in Daddi et al. (2007a), who 
assumed a linear correlation between fl-z color a nd E{B-V) 
derived from their SED fitting result ( IDaddi et al.ll2004) to es- 
timate extinction-corrected UV luminosities. In fact, if we 
apply the E(B-V) values estimated from the B-z color, 
the extinction corrected SFR(Ha) becomes smaller than the 
E{B-V) values estimated from our stellar population anal- 
ysis (Figure [TSll. Especially when the equal-extinction case 
is applied, the extinctio n corrected SFR(Ha)s are similar to 
the IDaddi eFal] (l2007al) relation. Although our E(B-V) val- 
ues correlat e with B-z color, th e correlation has larger scatter 
than that of IDaddi et all (l200l (Figure [IS]). The high-SSFR 
galaxies have larger E(B-V) values than the galaxies with 
similar fl-z color. 

The [Paddi et al. | (l2007al) relation is also supported by 
iPannella et al.l ( 120091) . who shows the averaged properties of 
extinction from the radio-to-UV luminosity ratio using a large 
numbers (^ 12000) of sBzK galaxies. They also find a strong 
correlation between the stellar mass and the amount of ex- 
tinction, and thus the correlation accentuates the Dad di et alj 
(|2()07a) relation. In Figure [17] the amount of extinction of 
our galaxies is shown as a function of the stellar mass. Our 
results show a larger scatter than the errors of the correla- 
tion by IPannella et alJ (l2009h . Although the extinction values 
of the low-SSFR gala xies corre late with the stellar masses, 
consistent with that of IPannella et alj (l2009l) . the high-SSFR 
galaxies have higher extinctions than the low-SSFR galax- 
ies with similar stellar masses. Because our emission line 
galaxies would be biased towards Ha bright systems, g alaxies 
t hat de viate from the average relation shown by Pannell a et alJ 
(120091) can be found, which probably results in the larger scat- 
ten 

IDaddi etal] (l2007al) also shows the same correlation us- 
ing a SFR indicator unrelated to extinction. They estimated 
the SFR from the sum of the extinction uncorrected UV and 
mid-IR luminosities. On the other hand, Cap uti et aP (l2006al) 
reported that there exist significant fraction of lower mass 
(- IO'^Mq) galaxies with high SSFRs ( > 10"^yr"') at z ^ 2 
using SFR inferred from the 24 /im flux. IDaddi etal] (l2007ah 
argues that the correlation is not accurately recovered with 
SFRs inferred only from mid-IR, due to a substantial pres- 
ence of MlPS-excess galaxies, which are excluded from their 
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Fig. 14. — Same as Figures fTol but plotte d with the results from literature. The black solid line shows the correlation between SFR and stellar mass of sBzK 
galaxies iDaddi et al. 2007a). The results of lErb et aH <2006bl) and lHavashi et alj <2009l) are shown by cyan circles and blue triangles, respectively. White dots 
denote MlPS-excess galaxies. 
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Fig. 15. — Same as Figures fTol but the amout of extinction is estimated from B—z color, as lDaddi et alj t2007af) . 



sample. A MlPS-excess galaxy could harbor a heavily ob- 
scured AGN (Daddi et al. 2007b). However, our results on 
the SFR inferred from Ha do not recover the correlation 
even if MlPS-excess galaxies are excluded from the sam- 
ple, ^^e_calculate_the MIPS excesses of our sample follow- 
ing lDaddi et alj (l2007ah . and found that only a fraction of the 
high-SSFR galaxies show MlPS-excess. In Figure[Tl MlPS- 
excess galaxies are denoted by white dots, which found to 
be more typically low-SSFR galaxies rather than high-SSFR 
galaxies. 

Previous studies on the relation between SFR and stellar 
mass using the SFR derived from the Ha luminosity show a 
similar trend to ours. The SFRs and stellar masses of BM/BX 
galaxies (Erb et al. 2006b) and sBzK galaxies (Hayashi et al. 
I2009h are also plotted on Figure [14] We convert the SFR 



and stellar mass of these studies to those calibrated with a 
Chabrier (2003) IMF In the literature, they applied differ- 
ent extinction corrections to estimate the Ha luminosities 
from that used to estimate the stellar extinction. The equal- 
extinction case is applied to the BM/BX galaxies, while the 
ICalzettil (1200 Ih relation to the sBzK galaxies. For compari- 
son, we applied the same extinction correction to the samples 
in the literature as that for the present sample in each panel. 
There are no correlations between SFR and stellar mass of the 
samples in the both panels. However, som e high-SSFR ga lax- 
ies with SSFR> 10"^ yr"' appear when the ICalzettil (120011) re- 
lation is applied. 

In the studies for BM/BX and sBzK galaxies, validity 
of the extinction correction is confirmed with the consis- 
tency between SFR(Ha) and SFR(UV) dErb et alj l2006bl: 
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Fig. 16. — Comparison of £(B — V) w ith B-z color for the present sam- 
ple. Symbols are same as in Figure 1121 The Dashed line shows the linear 
correlation between B—z color and E(B — V) bv lDaddi et alj i2004h . 



iHavashi et alj|2009l) . We also confirmed the consistency be- 
tween them in the previous section. The extinction corrected 
SFR(Hq!) values assuming the equal-extinction case are con- 
sistent only for lower SFR g alaxies (SFR< lOOMtnyr"'). 
SFRs calculated assuming the ICalzettil (120011) relation are 
broadly consistent, but the SFR(Ha)s are systematically 
larger than the SFR(UV)s by 0.3 dex. These are almost con- 
sistent with the previous SFR(Ha) works. Although Erb et al] 
(|2006b) applied only the equal-extinction case, galaxies with 
SFR( UV) > 100M(Ti yr"' are rai-e in their s ample. O n the other 
hand. IHavashi et alJ (2009). who used the ICalzettil (1200 lb re- 
lation, also showed that their SFR(Ha) values are systemati- 
cally larger than SFR(UV) by a factor of 3. 

Recently, the relation between Balmer decrement and stel- 
lar continuum attenuation h as been investigated for galaxies 
at low-to-medium redshift. ICowie & Barged ( 120081) studied 
spectroscopically confirmed A'-selected sample at 0.05 < z < 
0.475. Their extinctions from the Balmer decrement and the 
continuum extinctio ns are eq ual on average. On the other 
hand, Moustakas & Kennicutd(i2006 ) studied nearby galaxies, 
including starburst galaxies such as infrared-luminous galax- 
ies. Comparison between their extinctions shows large scatter, 
but the ave rage stellar-to -nebular extinction ratio is consistent 
with the Calzetti (2001) relation. The dust properties of the 
samples in these two studies are quite different. For exam- 
ple, the median reddening of iMoustakas & Kennicutj (120061) 
sa mple is E{B-V) ~C^24, while about half of the galaxies 
of ICowie & Bargerl ( l2008h have E(B-V) < 0.01. The values 
of reddening of the high-z samples discussed above and our 
samples are similar to that of Moustakas & Kennicutt (2006), 
but vary with the sample selection. The median values of 
reddening are E(B-V) - 0.15, 0.38 for BM/BX galaxies by 
lErb et al. (2006b) and sBzK galaxies by Havashi et al. (2009), 
respectively. This can imply that the iCalzettii (i2001i) relation 
gives better estimation of the value of reddening for ionized 
gas in the active star-forming galaxies at high redshift. 

5.2. Metallicity 

The gas-phase metallicity reflects the past star-formation 
history of a galaxy since the heavy elements are produced by 



Fig. 17. — Comparison of stella r ma ss with E(B — V) for the present sam- 
ple. Symbols are same as in Figure [T2l 



nuclear fusion in stars and fed back into gas. Therefore, it is 
expected that the gas-phase metallicity of galaxies in the high- 
z universe is lower than that in the local Universe. The metal- 
licities were found to be lower in high-z galaxies than those in 
local gala xies (e.g., Shaplev et al. 2004; Savaglio et al. 2005; 
Erb et all l20 06a: Maiolino et al. 2008; Havashi et al. 2009; 
Perez-Montero et al.,i2009; .Lamareille et al..i2009i) . We mea- 
sured the gas-phase metalUcity of our samples to investigate 
their star-formation histories. 

The gas-phase metallicity is estimated from the 
[N II] A6583-Ha line ratio (N2 index), defined by 



N2 = log{[Nn]A6583/Ha}, 



(9) 



which correlates with oxygen abundance 

( Storchi-Ber gmann et al. 1994). The relation between 
the N2 index and gas-phase oxygen abundance is given by 

12 + log(0/H) = 8.90 + 0.57 xA^2 (10) 

dPettini & Page! '2004). [Nil] and Ha are so close to each 
other in wavelength that the ratio is not affected by reddening. 
Therefore, the N2 index has been used to estimate the oxy- 
gen abu ndances of distant galaxies (e.g.. lShaplev et al]|2004t 
lErb et a l. 2006a; Hayashi et al. 2009). We note that the N2 in- 
dex sa turates at - IZq (12 + log(0/H) = 8.66; lAsplundet al.l 
I2OOI . 

In order to gain a sufficiently high S/N to reliably measure 
the line flux ratio, we stacked the [Nll]-Ha emission lines. 
We divided the present sample into three groups based on 
stellar masses at 9 x lO'^M© and 2.5 x lO'^M©. In addi- 
tion, we also divide the sample into two groups based on a 
SSFR of 10"** yr"'. The stacking, measurement of flux, and 
estimates of er ror are performed in the same manner as de- 
scribed in ^3.31 The results are plotted in Figure[T8]as a func- 
tion of the median stellar mass of each group. For compari- 
son, - 39500 local objects from SDSS DR-Cl local Lyman 
Break Analog (LBA) galaxies dOverzieret al.! 2010). B M/BX 
galax ies (.Erb et al...2006a) . and sBzK galaxies (.Havashi et alJ 
120091) are also plotted with the same calibration. 

'"^ Ihttp: // www.mpa-garching.mpg.cie/SDSS/DR4 /I 
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The mass-metallicity properties of our high-SSFR galax- 
ies are possibly different from that of rest-UV selected galax- 
ies. The mass-metallicity relation of the two largest mass 
bins is consistent with that of BM/BX galaxies described by 
lErb et al] (12006a'), while the smallest mass bin has highe r 
metalUcity than the similar mass bins of Er b et al.l (I2006ah : 
our sample is closer to the local galaxies. The discrepancy in 
the smallest mass bin could be due to the galaxies with high 
SSFR. The high-SSFR group has higher metallicity, close to 
that of the local galaxies with similar stellar mass. As shown 
in Figure[TO] high-SSFR galaxies dominate our sa mple in the 
smallest mass bin, while most of the galaxies in lErb et al.) 
(I2006ah have SSFR smaller than 10"** yr"'. 

The /T-band selected samples would provide many galax- 
ies with high SSFR and metallicity. The existence of higher 
metallicity galaxies in /ST-selected sample is also supported by 
the observations of sBzK galaxies by iHavashi et al.l (120091) . 
as shown in Figure [18] On the other hand, the local LBAs, 
which are local analogs of high-z UV-selected galaxies, tend 
to have a lower ma ss-metallicity relat ion. This would suggest 
that the sample by lErb et al.] (I2006al) is biased toward lower 
metallicities due to the rest-UV selection. The results indi- 
cate that in order to reveal the unbiased cosmic evolution of 
mass-metallicity relation, we need to evaluate metallicities of 
a larger sample of /T-band selected (i.e., mass-selected) galax- 
ies. 

It is surprising that the high-SSFR galaxies have already ac- 
quired metallicities as high as local galaxies at z '-^ 2. The age 
of the stellar population of the high-SSFR galaxies is young, 
so that the metal enrichment must have occurred in a short 
timescale. The N2 index is known to be sensitive to a hard 
ionizing radiation field from an AGN or stro ng shock exci- 
tation and ionization parameter (iKewlev & D opita 2002). In 
fact, our stacked [Olll]A5007/H;9 and [Nll]A6583/Ha line 
ratios show evidence for a harder ionizing spectrum. If the 
metallicity is really high, the result is in conflict with the ob- 
servations of local galaxies. These galaxies can not keep their 
stellar masses and metallicities after the observed epoch, be- 
cause the studies for local galaxies demonstrate that galax- 
ies in a similar mass range formed a s ignificant fraction o f 
their stars in the low redshift universe (iHeavens et al.ll2(J(J4l) . 
suggesting that these local galaxies must have increased their 
metallicity recently. It is possible that such high-SSFR galax- 
ies evolved into larger stellar mass galaxies with low star- 
forming efficiency in the local universe, depleting the gas 
and/or merging with higher mass galaxies. 

5.3. Contribution to the Cosmic SFR Density 

We estimate the contribution of the observed galaxies to 
the cosmic SFR density at z '-^ 2 with the SFR(Ha)s, divid- 
ing the sample into two groups with SSFR > 10~^ yr~' and 
< 10"* ^ yr"' for the extinction corrected SFRs with the Calzetti 
(120011) relation. Since we observed only a limited number of 
sBzK galaxies, we need to assume that the SFR and redshift 
distributions of the other galaxies are the same as the observed 
galaxies; this is not unreasonable since we randomly selected 
the present sample from the sBzK-MIPS and sBzK-non-MIPS 
population. First, we further divide the sample based on the 
MIPS detection in the public catalog. Second, we calculated 
the total SFR in a certain group; we multiply the total SFR 
in the group by the ratio of the number of the entire galax- 
ies in the group to that of the number of galaxies observed. 
In the calculation of total SFR, we assumed that the SFR of 
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Fig. 18. — Oxygen abundance derived from the N2 index for the stacked 
[Nn]-Ha spectra, which are grouped by (1) stellar mass at 9 X IO'Mq and 
2.5 X lO'^M© (black hexagons) and (2) SSFR at 10"* yr"' (black diamond 
and open square), as a function of the median stellar mass of each group. Hor- 
izontal en'or bars indicate the 64% quantiles of stellar masses in the respective 
bins, while vertical bars indicate the error of the N2 index measurements de- 
rived from bootstrap sampUng. The results of the stacked spectra of BM/BX 
galaxies by Erb et al. (2006a) are plotted with cyan circles. The sBzK galax- 
ies by Havashi et al. 1 2009) and the local Lyman Break Analog (LBA) galax- 
ies by Overzier et al. (2010) are plotted with blue downward triangles and red 
upward triangles, respectively. Gray contours show the plots of ~ 395000 
objects from SDSS DR4 with the same calibration of oxygen abundance. 
The horizontal dashed Hne shows solar metallicity (12 + log(0/H) = 8.66; 
Asplund et al. 2004), where the N2 index calibration saturates. 

the galaxies without Ha detection is negligible; we set those 
SFRs to zero. Then the total SFRs of galaxies with and with- 
out MIPS detection are summed. Third, we divide the to- 
tal SFR by th e ex pected detection rate of Ha emission line 
(87.7%; see ^4. lb to recover the Ha emission lines not de- 
tected due to strong sky emission lines or atmospheric absorp- 
tion bands. Finally, we divide the total SFRs by comoving 
volume of the survey region for 1.4 < z < 2.5. 
The resultant total cosmic SFR densities of the sBzK galax- 



ies are 0.136!;j:!gg 
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M(7^ yr~' Mpc for extinction 
correction calculated using the ICalzettil (1200 Ih relation and 
the equal-extinction case, respectively. As we see in the pre- 
vious discussions, the relation between extinction estimated 
from the stellar continuum and the emission lines probably 
depend on the intrinsic luminosities, and thus the SFRs. To 
exa mine thi s effect , when we apply extinction correction with 
the ICalzettil (120011) relation for the sBzK-MIPS galaxies and 
equal extinction for the sBzK-non-MIPS galaxies (differential 
correction), the result is 0.089;!;Qo*2r ^^ ^i^Y '^^^e of the extinc- 
tion correction, high-SSFR galaxies have a large contribution 
to the SFR density. Although these galaxies account for only 
r-^ 20% of all of the galaxies, they contribute to ^ 75% of 
the SF R density for the extinction correction with the ICalzettil 
(120011) relation and the differential extinction. Even for the 
equal-extinction case, they contribute to no less than ^ 50 % 
of the SFR density. 

We compare our SFR density at z ^ 2 with previous works 
(e.g.. lHopkinsl l2004tlHopkins & Beaco mll2006, and t he refer- 
ences therein; |Wang et al. 2006; Caput iet alj 120071) . If we 
use the same calibration of IMF as that of these studies, 
the SFR densities are 0.203, 0.044, O.134M0yr-'Mpc"^ for 



MOIRCS Deep Survey VI 



21 



the ICalzettil (1200 Ih relation, the equal-extinction case, and 
th e differential cor recti on, respectively. The values compiled 
bv lHopkinsI (l2004h and lHopkins &Beacoml (12006 ') are 0.05- 
0.38M(7) yr~'Mpc , wh ich are consistent with both of those 
with the ICalzettil ( 12001!) relation and differential extinction. 
Recent studies on the SFR density at z ~ 2 using IR luminosi- 
ties (850 /Ltm; IWang et al.1 120061: 24 /xm; ICaputi et al.1 l2007h 
show sHghtly lower SFR densities (0.08-0.09 Mq yr"' Mpc"^), 
which are still consistent with the differential extinction cor- 
rection. This would also suggest that the relation between 
the dust properties of the stellar continuum and the emission 
lines are different depending on the intrinsic SFR. The SFR 
density using the equal-extinction case is lower than that of 
any study. The lower SFR density implies that applying the 
equal-extinction case to all galaxies underestimates SFR(Ha). 

6. SUMMARY AND CONCLUSIONS 

We performed near-infrared spectroscopic observations for 
37 BzK color selected star-forming galaxies with MOIRCS 
on the Subaru Telescope. The sample galaxies are picked 
from the deep /Tj-selected catalog of MODS. We select 18 of 
our sample from the publicly available 24 //m source catalog 
of MIPS on board the Spitzer Space Telescope. In order to 
avoid contamination by AGNs, we exclude the BzKs detected 
in the Chandra Deep Field-North 2 Ms Point-Source Cata- 
logs. All of the sample galaxies have strong 1.6 /im bumps 
in their SEDs, so that their optical emission is thought to be 
dominated by the stellar continuum. The Ha emission lines 
are detected from 23 galaxies, of which the median redshift 
is z = 2.12. In the diagnostic diagram with [Olll]A5007/H/3 
and [N ll] A6583/Ha, the line ratios of the stacked spectra are 
above that found for local star-forming galaxies but closer to 
the expectations from star-formation activity rather than AGN 
activity. 

We derived the SFRs from the extinction-corrected Ha lu- 
minosities. The extinction correction for the Ha luminosity is 
estimated from the SED fitting of multi-band photometric data 
with stellar population model s of exponentially declining SFR 
and the Calzetti et al] (120001) extinction curve. To estimate 
the extinction correction for the ionized gas, we assumed the 
ICalzetti (2001) relation (£■(£- y)steiiar = 0.44 x E{B-V).^^^ as 
well as the equal-extinction case £(£- V)steiiar = EiB-V)^^^,^, 
and compared the extinction corrected SFRs. Although the 
Balmer de cremen t of the stacked emission lines is consistent 
with the C alzettil (|2001) relation, a fraction of the H/3 de- 
tected galaxies are consistent with the equal-extinction case. 
We further investigated these galaxies by comparing the ex- 
tinction corrected Ha luminosities with other SFR indicators. 
The comparison with the UV and MIPS 24 /xm flux shows 
that the SFRs of lower SFR galaxies (< lOOMgyr"') are in 
good agreement for the equal-extinc tion case , whil e those of 
the higher SFR galaxies agree for the lCalzettil (120011) relation. 
This implies that the relation between the dust properties of 
the stellar continuum and emission lines are different depend- 
ing on the intrinsic SFR (differential extinction correction). If 
the Calzetti relation is assumed, a fraction of the high-SFR 
galaxies should be detected in radio, (sub-)mm, and X-ray 
catalogs, but none of them except one AzTEC source are de- 
tected. We could not explain the reason for this discrepancy 
consistently. However, most of the discrepancy could be ex- 
plained by the short timescale of star-formation activity, be- 
cause the younger ages of the stellar population are derived 
from our SED analysis. We note that other mechanisms (e.g.. 



top-heavy IMF, lower dust temperature, etc.) could also ex- 
plain a part of the discrepancy. 

We compared the extinction-corrected Ha luminosities 
with the stellar masses obtained from the SED fitting. When 
we apply the Calzetti (2001) relation, the comparison shows 
no correlation between SFR and stellar mass. A fraction of the 
galaxies with stellar mass smaller than ~ lO'^Mo show SFRs 
higher than ~ lOOMQyr"^ The SSFRs of these galaxies are 
higher than ~ 10"'^ yr"'. Even if equal extinction is applied, 
their SSFRs are stifl high (SSFR > 3.5 x 10"\ r"'). This is 
inconsistent with the correlation demonstrated by lDaddi et al.l 
(2007a). Our value of extinction correction show s large scat- 
ter against B-z color and stellar mass, with which' Paddi et al.l 
(2004) and Pannella et al. (2009) show strong correlations. 
From the best-fit parameters of the SED fitting for these high- 
SSFR galaxies, we found that the average age of the stellar 
population is younger than 1 00 Myr, which are consistent with 
SFR age (SSFR"') with the lCalzettil (120011) relation. We also 
calculate the contribution of the present Ha emission to the 
cosmic SFR density at z '^ 2. Our results are consistent with 
previous works, if we apply the Calzetti. (.2001 ) relation or dif- 
ferential extinction correction. The high-SSFR galaxies sig- 
nificantly contribute (50-75%) to the cosmic SFR density. 

The metallicity of the high-SSFR galaxies, which is esti- 
mated from the N2 index, is larger than that expected from 
the mass-me t allicity relation of UV-selected galaxies at z ^ 2 
bv lErb et al.l (l2006al) . Although metallicities estimated from 
the N2 index is known to have significant uncertainties due to 
the harder ionizing radiation field and ionization parameter, 
the high metallicity comparable to galaxies in the same mass 
range is inconsistent with the known star-formation history of 
such galaxies in the local universe, which formed a significant 
fraction of their stars in the low redshift universe. The high- 
SSFR galaxies, therefore, possibly evolve into larger stellar 
mass galaxies with low star-forming efficiency, depleting the 
gas. 
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